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Abstract 

'To  assist  in  the  evaluation  of  proposed  structural  inspection  programs  for 
commercial  jet  transport  aircraft,  a logic  was  developed  to  simulate  structural 
defects,  failures,  and  Inspections.  This  logic  was  incorporated  In  a computer 
program  entitled  Structural  Area  Inspection  Frequency  Evaluation  (SAIFE) . With 
objective  of  quantifying  the  evaluation  process  currently  used  to  establish  and 
modify  inspection  intervals,  SAIFE  accounts  for  the  following  factors:  (1)  aircraft 
design  analysis;  (2)  fatigue  testing;  (3)  production,  service,  and  corrosion  de- 
fects; (4)  probability  of  crack  or  corrosion  detection;  and  (5)  aircraft  modifica- 
tion economics.  As  a five-volume  document,  this  report  covers  the  Initial  contract 
effort  plus  a subsequent  parametric  analysis  as  follows:  Volume  I (entitled  Execu- 
tive Summary)  presents  the  SAIFE  logic  and  documents  the  methodology  for  the  de- 
cision-making processes  in  the  simulation  logic.  Volume  II  details  the  SAIFE  sim- 
ulation logic,  presents  the  background  data  for  the  analytical  functions  and  de- 
cision-making processes,  and  includes  data  for  a typical  simulation  problem.  Vol- 
ume III  (entitled  Demonstration  Input,  Inspection  Survey,  and  MRR  Data)  presents 
data  tabulations  derived  from  historical  trends  and  design  Input  data  for  a SAIFE 
demonstration  problem.  As  the  user's  manual  for  the  SAIFE  computer  program,  Vol- 
ume IV, (entitled  Software  Documentation  and  User's  Manual)  contains  detailed  com- 
puter logic  flow  diagrams  and  a complete  listing  of  the  program  which  Is  written  In 
SIMSCRIPT  II. 5.  Volume  V (entitled  Results  of  Model  Demonstration)  presents  the 
results  of  the  program  application  to  a hypothetical  aircraft  and  compares  these 

results  with  the  service  experience  of  operational  aircraft. 
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TEMPESATURc 


PREFACE 


Technology  Incorporated  prepared  this  second  volume  of 
a five-volume  report  to  document  the  simulation  logic  for 
Structural  Area  Inspection  Frequency  Evaluation  in  accord- 
ance with  Article  II,  paragraph  B of  Contract  DOT-FA74WA- 
3493.  (Volume  II  along  with  Volume  I completes  the  require- 
ments of  Phases  I and  II  of  the  contract.)  The  effort  is 
sponsored  by  the  Aircraft  Safety  and  Noise  Abatement  Divi- 
sion, Systems  Research  and  Development  Service  of  the  Fed- 
eral Aviation  Administration. 

The  principal  Technology  Incorporated  personnel  engaged 
on  this  program  were  Mr.  Carter  J,  Dinkeloo,  project  engi- 
neer, who  served  as  principal  investigator;  Mr.  Martin  S. 
Moran,  research  engineer,  who  developed  the  model  for  the 
SAIFE  computer  program;  and  Mr.  Ronald  I.  Rockafellow,  pro- 
gram manager. 

The  contract  monitors  for  the  FAA  were  Messrs.  Herbert 
Spicer  and  Charles  Troha  of  the  Aircraft  Safety  and  Noise 
Abatement  Division.  The  technical  monitor  was  Mr.  Arnold  li. 
Anderjaska  of  the  Flight  Standards  Division. 
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INTRODUCTION 


I . 


It  is  the  mutual  goal  of  the  FAA , airframe  manufactur- 
ers, and  air  carriers  to  constantly  improve  the  structural 
integrity  and  inspection  efficiency  of  civil  aircraft.  The 
good  safety  record  of  U.S.  air  carriers  indicates  that  the 
current  process  of  establishing  and  modifying  structural 
inspection  programs  has  boon  successful.  However,  with  the 
increasing  size  and  complexity  of  second-  and  third-genera- 
tion transport  aircraft,  there  is  a need  to  quantify  more 
precisely  the  current  subjective  evaluation  process  which 
relies  heavily  on  reliability  analyses  of  the  new  design  and 
on  operational  experience  of  similar  aircraft. 

because  of  the  extreme  complexity  of  the  evaluation 
process,  a computer  simulation  of  all  critical  aircraft 
service  life  aspects  was  judged  the  most  rational  means  for 
quantifying  the  process  more  exactly,  As  a five-volume 
document,  this  report  documents  the  resultant  Structural 
Area  Inspection  Frequency  livn  1 uat  Lon  (SAIITil  simulation 
logic.  S A 1 1 ■ 1 i accounts  for  the  following  factors:  (.1)  air- 
craft design  analysis;  (2)  component  and  full-scale  fatigue 
testing;  (.3)  production,  service,  and  corrosion  defects; 

(4)  probability  of  crack  or  corrosion  detection;  and  (5)  air- 
craft modification  economics.  It  treats  these  factors  in  a 
logical  sequence  that  realistically  represents  the  procedure 
currently  used  to  establish  and  modify  inspection  intervals. 
SAIF ti  is  designed  to  provide  a repeatable  method  for  evalu- 
ating proposed  inspection  programs.  However,  it  is  not 
intended  to  supplant  the  Maintenance  Review  Hoard  or  the  air 
carrier  use  of  the  Standard  Operations  Specification  - 
Aircraft  Maintenance . 

In  addition  to  presenting  the  SAIF li  logic  applicable  to 
initial  demonstration,  this  report  documents  the  research 
conducted  to  establish  the  quantitative  functions  required 
for  decision  logic  in  the  simulation.  Some  of  the  documen- 
tation for  these  functions,  such  as  fatigue  life  scatter, 
are  taken  from  work  conducted  in  other  studies.  Other 
functions,  such  as  the  probability  of  defect  detection,  are 
the  result  of  work  conducted  as  part  of  this  contract. 
Whatever  the  source,  all  analytical  information  is  refer- 
enced throughout  the  report.  The  logic  applicable  to  the 
parametric  study  is  given  in  Volume  IV,  Hook  2 of  this 
report . 

Figure  1 illustrates  the  data  sources  ami  analytical 
functions  that  are  integrated  into  the  SAIFH  logic.  As 
Volume  11,  this  volume  presents  the  detailed  simulation 
logic  incorporated  in  SAlFli  and  all  the  background  data 
required  for  the  analytical  functions  and  decision-making 
processes.  It  also  includes  much  of  the  data  required  for  a 
typical  simulation  program. 


II.  DISCUSSION  OF  DETAILED  SIMULATION  LOGIC 


The  eight  blocks  in  Figure  2 represent  the  major  as- 
pects of  the  SAIFE  simulation  logic  which  was  originally 


Figure  2.  Flow  Diagram  Showing  Major  Aspects  of  SAIFE  Logic 


3 


developed  by  Anderjaska  in  Reference  1.  block  1.0  accepts 
input  data  for  the  aircraft  fleet  and  for  each  structural 
element  in  these  aircraft.  After  determining  whether  ele- 
ment modifications  are  required  because  of  the  fatigue  test 
results  in  Block  2.0,  Block  1,(1  assigns  a fatigue  life  to 
each  element  in  each  aircraft.  Block  3.0  determines  whether 
production,  service,  or  corrosion  defects  will  occur;  if  it 
is  determined  that  such  defects  will  occur,  Block  3.0 
predicts  the  times  when  they  will  occur.  After  comparing 
the  flight  loads  with  the  strength  of  each  element,  Block 
4,0  predicts  the  time  to  structural  failure  for  each  air- 
craft. Block  5.0  conducts  the  periodic  inspections  of  each 
element.  If  a structural  failure  has  occurred.  Block  5.0 
deletes  the  corresponding  aircraft  from  the  fleet.  However, 
if  an  element  has  a defect  that  is  detected.  Block  0,0 
repairs  the  element.  Depending  on  the  magnitude  of  the 
detected  defects,  special  inspections  and  increased  inspec- 
tion frequencies  may  be  called  for  in  Block  7.0  and  modifi- 
cations may  be  instituted  in  Block  8.0.  When  all  the  air- 
craft have  boon  deleted  from  the  fleet  either  through  re- 
tirement from  service  or  as  a result  of  a structural  failure, 
the  simulation  is  complete. 

As  detailed  in  the  following  sections,  each  block  is 
broken  down  into  several  subblocks  which  contain  the  func- 
tions and  decision  logic  required  for  the  simulation,  liach 
subblock  for  a particular  block  breakdown  is  represented  by 
the  unit,  number  Identifying  the  block  in  Figure  2 and  by  a" 
decimal  number  denoting  the  subblock  sequence;  for  example, 
the  subblock  represented  by  2.3  correlates  with  Block  2.0  In 
Figure  2 and  is  third  in  the  subblock  sequence.  By  design, 
the  subsection  numbering  in  this  section  correlates  with 
this  numbering  system. 
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Detailed  Description  of  Block  1.0,  Input  Data/Generate 
Fatigue  Lives  (see  Figure  3) 


Figure  3.  Detailed  Flow  Diagram  - Input  Data/Generate 
Fatigue  Lives 


1 • 1 Input  Aircraft  Fleet  Data 

This  block  introduces  the  following  input  data  which  is 
constant:  Tor  all  elements  throughout  the  simulation. 

(1)  Aircraft  type  (two-clement  array). 

(2)  The  size  of  fleet  to  he  simulated. 

(3)  Service  life  of  the  aircraft  (flight  hours). 

(4)  Time  ut  which  aircraft  production  begins  (flight 
hours  after  start  of  simulation). 

(5)  Initial  aircraft  production  (flight  hours  between 
aircraft). 

(6)  Second  aircraft  production  rate  (flight  hours 
between  aircraft). 

(7)  Time  at  which  second  production  rate  takes  effect 
(flight  hours  after  start  of  production). 

(8)  Time  at  which  fatigue  test  is  started  (flight 
hours  after  start  of  simulation). 

(9)  Fatigue  test  acceleration  factor. 

(10)  Corrosion  area  growth  rate  (square  Inches  per 
•flight  hour). 

(11)  Percentage  of  fatigue  life  at  which  inspection 
frequency  is  increased  because  of  a fatigue  test 
fa i lure . 

(12)  Factor  to  reduce  fatigue  life  when  corrosion  oc- 
curs in  a stress  concentration. 

(13)  Factor  to  reduce  fatigue  life  when  corrosion  oc- 
curs outside  n stress  concentration. 

(14)  Mean  (iig)  and  standard  deviation  (eg)  for  the  log- 
normal  distribution  of  the  ratio  of  the  actual 
average  fatigue  life  to  the  predicted  average 
fatigue  life. 

(15)  liquation  constants  (two  required)  for  gust  and 
maneuver  load  distribution. 

(16)  Initial  inspection  intervals  (four  required) 
(flight  hours). 

(17)  Cost  of  each  level  of  inspection  (four  required) 
(dollars) . 


(3 


(18)  D-level  sampling  percentage. 

(19)  A "YliS"  or  "NO"  switch  used  to  implement  the  long 

list  output  option. 

(20)  Percentage  of  critical  crack  length  at  which  an 
internal  crack  becomes  external. 

A total  of  29  input  constants  are  required  for  this 
data  block. 

1 . 2 Input  Element  Data 

This  block  introduces  the  following  input  data  which  is 
unique  to  each  element  (therefore,  each  element  must  have  a 
complete  set  of  data  specific  to  the  given  element): 

(1)  Element  identification  (four-element  array). 

(2)  Predicted  average  fatigue  life  (flight  hours), 

(3)  Actual  average  fatigue  life  (flight  hours). 

(4)  Average  slow  crack  growth  rate  (inches  per  flight 
hour) . 

(5)  Average  fast  crack  growth  rate  (inches  per  flight 
hour) . 

(6)  Crack  length  to  structural  failure  (inches). 

(7)  Critical  crack  length  (inches). 

(8)  Pail -safe  crack  length  (inches). 

(9)  Probability  of  production  defect  (occurrences  per 
aircraft) . 

(10)  Corrosion  resistance  rating. 

(11)  Service  damage  occurrence  rate  (occurrences  per 
flight  hour) . 

(12)  Lead  time  to  implement  modifications  (equivalent 
hours) . 

(13)  Factor  by  which  inspection  intervals  arc  decreased 
because  of  fatigue  test  failure  (applied  to  "C" 
and  "D"  level  inspections  only). 

(14)  Factor  by  which  inspection  intervals  arc  decreased 
because  of  unfavorable  service  experience  (applies 
to  "C"  and  "D"  level  inspections  only). 
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(15)  Factor  by  which  inspection  intervals  are  increased 
because  of  favorable  service  experience  (applies 
to  "C"  and  "D"  level  inspections  only). 

(16)  Probability  of  fatigue  crack  initiating  internally, 

(17)  Probability  of  corrosion  initiating  internally. 

(18)  Lowest  internal  inspection  level. 

(.19)  Lowest  external  inspection  level. 

(20)  Repair  costs  for  each  level  of  inspection  (four 
constants  required)  (dollars). 

(21)  Decision  variable  indicating  whether  or  not  modi- 
fications are  to  be  fatigue-tested  (yes  or  no). 

(22)  First  modification  tooling  cost  (dollars). 

(23)  Additional  modification  tooling  cost  (dollars). 

(24)  First  modification  installation  cost  (dollars). 

(25)  Additional  modification  installation  cost  (dollars). 

(26)  Repair  cost  for  defect  found  during  special  In- 
spection (dollars). 

(27)  Probability  of  existing  corrosion  being  in  a 
stress  concentration. 

(28)  Initial  corrosion  occurrence  rate  (occurrences  per 
flight  hour). 

(29)  Second  corrosion  occurrence  rate  (occurrences  per 
flight  hour). 

(30)  Airframe  time  at  which  second  rate  takes  effect 
(flight  hours). 

A total  of  36  input  constants  nre  required  for  this 
data  block. 

bach  element  to  be  evaluated  by  the  simulation  is 
identified  by  three  groups  of  alpha  characters  and  by  one 
group  of  numeric  characters.  The  alpha  characters  define 
the  basic  element  type  and  the  general  location  on  the 
aircraft,  while  the  numeric  characters  define  the  specific 
location  of  the  element  by  Identifying  the  wing  or  fuselage 
station  number.  For  example,  an  element.  Identified  as  "FUS- 
MFR-TOP-400"  would  be  a frame  locuted  in  the  fuselage  crown 
with  the  attaching  structure  extending  from  station  390  to 
station  410. 
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The  simulation  is  designed  to  handle  as  many  individual 
elements  in  each  aircraft  as  is  necessary.  Accordingly,  the 
size  of  each  wing  or  fuselage  clement  depends  only  on  natur- 
ally occur i ing  design  points  such  as  rib  or  fruinu  spacing. 
Therefore,  tht  element  identified  at  fuselage  station  400 
includes  all  structure  and  attaching  parts  between  fuselage 
stations  390  and  410,  In  this  example  the  fuselage  frame 
element  would  also  include  all  the  attached  skin  as  shown  in 
Figure  4.  This  figure  also  shows  a typical  wing  stringer 
element  witli  attaching  structure.  The  basic  element  types 
and  the  number  of  individual  elements  in  each  basic  type  arc 
listed  in  Table  1.  This  table  Is  applicable  to  a typical 
narrow-body  aircraft  and  is  used  throughout  this  report  to 
analyze  service  history  data  from  narrow-body  aircraft.  The 
identification  system  is  the  same  as  that  used  to  process 
the  MRR/SDR  historical  data  (Volume  III).  It.  offers  a great 
deal  of  flexibility  in  laying  out  the  elements  on  any  par- 
ticular aircraft  and  permits  an  easy  comparison  of  the 
simulation  output  and  the  historical  information. 
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Figure  4.  Typical  Airframe  Elements  Used  in  the  SAIF1I 
Simulation 


9 


TABU;  1.  ELEMENT  Dli EINlTtON  AND  DISTRIBUTION  I:OR 
TYPICAL  NARROW-BODY  JUT  TRANSPORT 

No.  Elements  in 

Basic  Typo  Bach  Aircraft 


WNG  (wing)  : 


STR  (stringer)  - 

FWD  (forward) 

S3 

- 

CRN  (center) 

33 

AFT  (aft) 

33 

MB  (rib) 

FW1)  (forward) 

33 

- 

CRN  (center) 

33 

- 

AFT  (aft) 

33 

SPR  (spar) 

RW1)  (forward) 

3 3 

- 

CRN  (center) 

33 

* 

ART  (aft) 

33 

ACC  (access) 

RRM  (f runic) 

50 

WSC  (wing  center  section) ; 

SWB  (spanwise  beam)  - RWD  (forward) 

7 

- CRN  (center) 

7 

-•  AFT  (aft) 

7 

RIB  (rib) 

• FWD  (forward) 

7 

- CRN  (eentor) 

7 

■ AFT  (aft) 

7 

STR  (stringer) 

■ FWD  (forward) 

7 

■ CRN  (center) 

7 

- ART  (aft) 

7 

FUS  (fuselage): 

MFR  (main  frame) 

- TOP  (top) 

60 

■ SID  (side) 

60 

- BOT  (bottom) 

60 

STR  (stringer) 

- TOP  (top) 

60 

- SID  (side) 

60 

- BOT  (bottom) 

60 

FLR  (floor) 

- BUM  (beam) 

6 0 

KliL  (keel) 

- BUM  (beam) 

6 0 

WIN  (window) 

- RRM  (frame) 

SO 

DOR  (door) 

- RRM  (Frnmc) 

10 

PRS  (pressure) 

- WI.B  (web) 

60 

1 . 3 End  of  Data 

This  function  simply  monitors  the  input  data  for  the 
end  of  data  card.  When  the  appropriate  end  of  data  code  is 
encountered,  the  simulation  is  concluded. 
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1. 4 Predicted  Average  Fatigue  Life 

A program  input,  this  variable  is  the  average  fatigue 
life  analytically  predicted  for  each  clement  design. 

As  used  in  this  report,  fatigue  life  is  defined  as  the 
accumulated  operational  time  when  a crack  initiates.  SAIFF. 
uses  the  following  three  types  of  fatigue  lives  which  are 
used  throughout  this  report: 

(1)  Predicted  average  fatigue  life  - the  fleet  average 
fatigue  life  for  an  element ’ determined  from  the 
manufacturer's  design  analysis. 

(2)  Actual  average  fatigue  life  - the  fleet  average 
fatigue  life  for  an  element  determined  from  ex- 
tensive service  experience. 

(3)  Actual  element  fatigue  life  - the  fatigue  life  of 
an  element  on  an  individual  aircraft  in  the  fleet. 

The  primary  source  of  data  for  this  variable  is  the 
manufacturer's  design  analysis.  When  such  data  is  not 
available  directly  from  the  manufacturer,  it  can  be  obtained 
in  a less  detailed  format  from  a Maintenance  Review  Hoard 
report  or  a Fatigue  Integrity  Program  report.  It  is  also 
possible  to  calculate  the  fatigue  life  from  service  experi- 
ence by  using  the  method  presented  in  Reference  2.  Tf  none 
of  these  sources  arc  available,  the  average  fatigue  life 
must  be  approximated  on  the  basis  of  the  design  service  life 
for  the  current  aircraft  or  from  a previous  aircraft  of  the 
same  manufacturer. 

1 • 5 Actual  Average  Fatigue  fife 

Since  fatigue  phenomena  are  not  completely  defined,  the 
fatigue  life  prediction  analysis  should  be  performed  statis- 
tically. Since  in  practice  a statistical  approach  is  not 
used,  the  actual  fatigue  life  of  n structure  of  a given 
design  will  usually  differ  from  that  analytically  predicted. 
The  probability  of  the  actual  life  being  greater  or  less 
than  that  predicted  was  studied  by  the  Royal  Aircraft  li  s - 
tablishment  (Reference  3).  In  Reference  3,  fatigue  lives 
based  on  full-scale  structural  fatigue  tests  on  the  wings  of 
British  military  and  civil  aircraft  were  compared  with  those 
based  on  average  fatigue  performance  In  laboratory  tests  of 
typical  aircraft  joints.  In  each  type  of  fatigue  life 
derivation,  the  estimated  life  was  based  on  average  fatigue 
performance  and,  as  necessary,  on  Miner's  linear  cumulative 
damage  method.  As  shown  in  Figure  5,  the  results  of  this 
comparison  indicate  that  the  actual  life  may  frequently  he 
overestimated  if  the  calculated  fatigue  performance  Is  not 
statistically  interpreted.  To  determine  the  distribution 
shown  in  Figure  5,  a computer  program  was  used  to  fit  the 
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following  statistical  ( 1 og -no  rma  I ) distribution  to  the  data: 


* 


„ _ Actual  Life 
R PTeTTFcOTTiTe 


(1) 


M ^ ■ 0.841 


■ 0.695 

The  parameters  pg  (distribution  mean)  and  og  (distri- 
bution standard  deviation)  are  input  variables.  These  pa- 
rameters enable  SAIFK  to  account  for  improvements  in  fatigue 
analysis  techniques.  An  example  of  the  relationship  result- 
ing from  improved  analysis  techniques  is  shown  in  Figure  5, 
where  up  ■ 1.000  and  eg  ■ 0.695. 


Figure  5.  Comparison  of  Predicted  and  Actual  Fatigue  Life 
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With  the  application  of  Monte  Carlo  computer  tech- 
niques, a log-normally  distributed  correction  factor  for  the 
predicted  average  fatigue  life  (block  1.4,  figure  .1)  is 
generated  for  each  element  design.  This  correction  factor 
(R)  yields  the  actual  average  fatigue  life  (block  1.4, 

Figure  3)  and  is  calculated  for  each  element  . The  same 
statistical  model  is  used  for  all  the  elements  making  up  an 
aircraft  model. 

1 . 6 Production  Modification  Pending 

When  each  aircraft  enters  servicb,  this  function  deter- 
mines whether  previous  simulation  logic  has  Instituted  a 
modification  because  of  operational  experience  on  aircraft 
already  in  service  or  because  of  a previous  fatigue  test 
failure.  If  a modification  is  pending,  the  logic  proceeds 
to  Block  1.7  to  determine  whether  the  modification  is  avail- 
able for  installation.  If  a modification  is  not  pending, 
the  logic  goes  to  the  routine  TMOD.  It  is  assumed  that  if 
previous  unfavorable  operational  experience  has  caused  the 
development  of  an  element  modification,  a subsequent  fatigue 
test  failure  will  not  cause  further  modification  of  that 
element. 

1.7  Modification  Installed 


This  function  indicates  whether  or  not  a proscribed 
structural  modification  has  been  installed  in  an  aircraft 
before  it  enters  service.  The  modification  has  been  in- 
stalled if 


X'l  * C2  :■  X 8 + fll  + Cl 2 (2) 


where  X4  3 production  number  of  aircraft 

C2  ■ production  rate  of  aircraft 

X8  ■ time  between  the  delivery  of  the  first  air- 
craft and  the  time  of  the  decision  to  develop 
a structural  modification  for  a given  element 

Cll  ■ lead  time  to  develop  a structural  modifica- 
tion for  a given  element  type 

C L 2 * load  time  from  the  development  to  the  produc- 
tion-line incorporation  of  the  structural 
modification  for  a given  element  type 

Items  C 2 , Cll,  and  C12  are  program  inputs.  Item  U2 
depends  on  the  aircraft  type  being  considered,  while  Cll  and 
Cl  2 depend  on  the  particular  element  being  considered.  Any 
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data  available  only  in  terms  of  calendar  days  must  be  con- 
verted to  the  equivalent  simulation  time  (flight  hours). 

Calendar  days  * a . 2 - equivalent  simulation  time  (3) 

This  conversion  is  based  oil  the  average  yearly  aircraft 
utilization  of  3000  flight  hours  per  365  days. 

The  modification  will  be  installed  on  aircraft  that  are 
already  in  service  when  the  aircraft  are  out  of  service  for 
repairs  or  for  normally  scheduled  overhaul  inspections  if  it 
has  been  determined  that  retrofit  modifications  are  economi- 
cally feasible  or  required  for  safety* of -flight  reasons. 

The  validity  of  this  function  depends  on  the  accuracy 
of  Cll  and  C12.  The  best  estimation  of  lead  times  should 
come  from  the  aircraft  manufacturer  or  from  the  air  carrier 
if  the  carrier  performs  the  modification.  The  lead  time  for 
a critical  item  that  precedes  the  modification  may  be  pro- 
vided by  a vendor  or  independent  supplier. 

1 . 8 Actual  11  Lenient  Fatigue  Life 

If  identical  fatigue  tests  are  performed  on  several 
nominally  identical  test  specimens,  the  resulting  fatigue 
lives  will  not  be  identical.  This  basic  fatigue  life  scat- 
ter is  a function  of  the  material  properties,  manufacturing 
quality,  and  process  variations  of  the  test  specimens. 

Also,  no  two  aircraft  within  the  same  fleet  experience 
identical  load  spectra.  This  load  environment  variation 
introduces  additional  fatigue  life  scatter  among  like  struc- 
tural elements  within  the  same  fleet. 

Clearly  then,  the  fatigue  life  of  aircraft  structures 
must,  be  treated  as  a stochastic  variable  whose  frequency 
distribution  reflects  both  the  basic  fatigue  scatter  and  the 
loud  environment  variation.  Because,  in  the  analysis  of 
aircraft  structures,  only  a small  sample  can  be  tested  and 
the  failure  of  even  a single  structure  may  be  catastrophic, 
the  expected  time  to  first  failure  is  a much  more  signifi- 
cant measure  than  the  mean  time  to  failure.  A method  for 
the  estimate  of  the  expected  time  to  first  failure  is  out- 
lined by  Freudenthal  (Reference  4).  in  his  paper,  Freuden- 
thal shows  that  the  cumulative  distribution  function  has  the 
form  of  the  two -parameter  Weibull: 

Fit)  « 1 - expf-  (t/(-))h]  (4) 


where  t = time  to  crack  initiation 
0 a characteristic  value 
b = shape  parameter 
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By  pooling  the  results  from  various  sources  (Figure  6), 
Freudenthal  makes  an  estimate  of  fatigue  life  scatter  based 
on  the  representation  of  test  data  by  the  log-normal  dis- 
tribution, He  concludes  that  a value  of  the  standard  devia- 
tion cr(logioN)  “ 0.15  - 0.20  is  representative  of  most 
results  in  the  long-life  range  (N  > 10&  cycles).  These 
values  are  consistent  with  the  work  done  by  Abelkis  (Ref- 
erence 5)  in  which  he  concludes  that  a value  of  o = 0.14 
would  describe  the  basic  fatigue  life  scatter  and  that  a 
value  of  a ■ 0.20  would  account  for  basic  scatter  and  the 
additional  scatter  introduced  by  load  environment  variation. 
Using  the  actual  average  fatigue  life'(t)  generated  by  Block 
1.5  and  assuming  a value  for  a,  the  values  of  0 and  b can  be 
determined  from  the  following  relationship: 


n/(2.505  a /%  ) 


t/ro  + 1/b) 


where  f ■ the  gamma  function 


• AUSTRALIAN  AIRCRAFT  WING  TESTS 
O COLUMBIA  UNIVERSITY  RANDOM  LOAD  TESTS 
'o  COLUMBIA  UNIVERSITY  CONSTANT  LOAD  TESTS 
A NLL  AMSTERDAM  PROGRAM  LOAD  TESTS 


I 

1 


25  3.0  35  AO  A5  t ^ 5,0  55  6,0  65  7,0 


Figure  6.  Relation  Between  6 - o(loglf)  N)  and  log  A 
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The  Monte  Carlo  technique  can  now  be  applied  to  generate  a 
two -parameter  Wei  bull  distribution  of  fatigue  lives  by 
replacing  [l-l:(t)|  in  liquation  H)  with  a random  number  for 
each  aircraft  in  the  .simulated  fleet; 

t « 0 [ In ( 1/RN) ] 1/b  (7) 


where  RN  * a uniformly  distributed  random  number. 

The  simulation  program  allows  any  structural  element  to 
have  as  many  as  three  fatigue  cracks  at  the  same  time.  The 
times  to  subsequent  crack  initiations  are  arrived  at  by 
truncating  the  distribution  at  the  previous  fatigue  life  and 
drawing  the  time  to  the  next  crack  initiation  from  the  new 
conditional  distribution.  Let  Fj(t)  be  the  cumulative 
distribution  function  of  time  to  crack  initiation.  Make  a 
random  draw  from  Fi(t)  which  yields  ti  (time  to  first  crack 
initiation').  Now  truncate  the  lower  end  of  the  distribution 
at  ti  and  make  a random  draw  from  the  new  conditional  dis- 
tribution 1; 2 ( t ) to  yield  ti  (time  to  second  crack  initia- 
tion) where 


Lift)  - r- 1 ( t ! ) 

I-  2 ( t ) ■ Hi  f t 1 1 > 1 1 ) ■ — — ( 8 ) 


and  P,(t|t  > t.)  ■ probability  of  crack  initiation  before 

time  t with  the  condition  that  crack 
initiation  occurs  after  time 

F,  (t)  ■ probability  of  crack  initiation  before 
time  t 


P 1 C 1 1 ) a probability  of  crack  initiation  before 
time  t;( 

1 - IT i ( 1 1 ) » probability  of  no  crack  in  i tiation  before 
time  t i 


Now  truncate  the  lower  end  of  Fjft)  at  tz  and  make  a 
random  draw  from  the  new  conditional  distribution  F3U)  to 
yield  tj  (time  to  third  crack  initiation)  where 


F3(t) 


1;2  ft: ! t > t2)  ■ 


r; 2 c t ) - r,ct2) 

““T^-TTr tp— 


(9) 


How  the  three  probability  density  functions  f 1 ( t ) f 
1 2 f t ) » and  f 3 ft)  are  related  can  be  determined  from  their 
definitions.'  From  probability  theory, 
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i. . -y...  • -r,r  n . r, 


lkaamfifeafaiiaii tauttiAvwHUA^vtt..  intaiitomnahPihiT -t  . ■i.vr  ■ > „ > ■ . in  .nj,.  u* 


(10) 


fx(t) 


dFj  (t) 

3t 


Likewise, 


f2(t) 


dP2(t)  d FjCt)  - Pl(tl) 

at  T-T1Tt1) 


or 

MO  , 

f2(t)  ■ r~-F',t't1y 

Similarly, 

r , MO  f ! CtD 

f3(t)  " rfjrrp' " irr^rpi  [t-p2 ct~TT 


(id 


(12) 


(13) 


If  we  let  Ki  ■ 1/ ( 1 ‘Fi (ti ) ) and  K2 
can  then  write  f - •*  - -•  * 
and  ftft)  can  be 


. and  K2  - 1/[1-Fi(ti)] fl-F2(t2) j , 

MO  " KifiCt),  f3(t)  - K 2 f 1 ( t J . Thus  f2(t) 
and  ft(t)  can  be  formed  from  fi(t)  by  truncating  the  lower 
end  at  ti  and  t2  and  by  multiplying  all  values  of  £i(t)  by 
the  constants  Kj  and  K2,  respectively.  This  process  is 


demonstrated  in  Figure  7. 


Figure  7.  Truncated  Fatigue  Life  Distribution  for  First, 

Second,  and  Third  Cracks  in  an  Individual  Element 


we 
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In  Figure  8,  the  cumulative  distribution  ot  time  to 
crack  Initiation  is  shown  for  each  of  the  three  cracks  on 
the  individual  elements.  Const  rue.  ted  from  the  results  of  a 
sample  simulation  problem,  those  distributions  show  a dis- 
tinct Weibull  form  w i t h prog  ress  i ve  1 >’4  higher  characteristic 
vaiues  10).  This  is  the  expected  result  from  the  method 
discussed  above. 


< 

J ,\ 


I’l.KiHT  T I Mis  (loon  lltHJUS) 

Figure  8.  Cumulative  Distribution  of  Fatigue  Lives  for 
First,  Second,  and  Third  Cracks 
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1 . 9 Aircraft  Enters  Service 

This  function  causes  aircraft  to  enter  service  at  a 
production  rate  prescribed  in  the  input.  The  program  logic 
has  the  facility  to  change  the  production  rate  once  at  some 
time  during  production.  The  time  of  this  change  and  the  new 
production  rate  are  also  prescribed  in  the  input.  Aircraft 
continue  to  enter  service  until  the  fleet  size  defined  in 
the  input  is  attained. 


DECISION  HADE  TO 
DEVELOP  MODIFICATION 


INSPECTION  FREQUENCY 
INCREASED  PENDING 
MODIFICATION 


MODIFICATION  TESTED 
(INPUT) 


ACTUAL  AVERAGE  MODIFIED 
FATIGUE  LIFE 


Figure  9,  Detailed  Flow  Diagram  ■ Develop  Modification 
Because  of  Fatigue  Test  Failure 

2 . 1 Decision  Made  to  Develop  Modification 

This  function  determines  whether  a modification  to  an 
element  is  required  so  that  the  aircraft  type  may  reach  its 
predicted  service  life.  The  decision  is  based  on  the  number 
of  equivalent  flight  hours  attained  during  fatigue  testing. 
For  airframe  elements,  a goal  of  two  times  the  service  life 
is  commonly  used.  If  this  goal  is  achieved,  the  fatigue 
test  may  be  discontinued  or  it  may  be  continued  to  determine 
what  additional  margin  of  safety  is  present. 


To  dotermino  whether  the  fatigue  test  goal  has  been 
achieved,  the  hours  of  testing  are  multiplied  as  follows 
by  a fatigue  test  acceleration  factor  to  arrive  at  the 

equivalent  flight  hours: 

Test  Hours  * Fatigue  Test  Acceleration  Factor 

■ Flight  Hours  (14) 

The  criteria  for  developing  u modification  are  then  as 
follows: 

(1)  If 

flight  hours  _>  2 * service  life  (15) 

a modification  is  not  developed. 


(2)  If 

flight  hours  < ?,  * service  life  (16) 

a modification  is  developed  and  it  is  in- 
stalled at  production  when  it  becomes 
available . 

Because  of  the  significantly  higher  cost  of  installing 
a modification  on  an  aircraft  already  In  service,  retrofit 
modifications  arc  not  Installed  unless  the  fatigue  test 
failure  occurred  in  loss  than  one  service  Life,  If  the  fa- 
tigue test  failure  occurred  In  less  than  one  service  life, 
the  modification  Is  required  for  safotv-of- f 1 ight  and  is 
installed  on  all  aircraft. 

2 . 2 Inspection  Frequency  Increased  Ponding  Modification 

When  it  has  been  determined  that  a modification  must 
be  installed,  there  is  a lead  time  required  to  design  and 
fabricate  the  modification  and  to  await  the  aircraft's  being 
scheduled  -for  an  out -o f- service  period.  During  this  lead 
time,  it  may  be  necessary  to  Increase  the  frequencies  of  the 
lowest  level  close  internal  and  close  external  inspections. 

The  decision  to  increase  the  Inspection  frequency  is  based 
on  the  assumption  that  the  fatigue  test  specimen  represent!: 
an  average  of  all  elements  and  that  a scatter  factor  is 
required  to  account  for  nil  the  elements  In  a typical  fa- 
tigue life  distribution.  Therefore,  when  the  flight  hours 
on  any  particular  aircraft  roach  some  percentage  of  the 
fatigue  test  failure  life,  either  that,  aircraft  must  be 
modified  or  the  inspection  frequency  must  ho  increased  until 
tho  modification  is  installed.  The  percentage  of  fatigue 
test  failure  life  at  which  the  inspection  frequency  is 
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increased  is  ari  input  parameter;  in  addition,  the  factor 
which  increases  the  inspection  frequency  for  both  the  close 
external  and  the  close  internal  inspection  is  also  an  input 
etrameter  since  r t ■ a required  change  depends  on  the  element 
! ru>  considered. 

- • Modification  Tested 

When  block  2.1  Indicates  that  a modification  is  re- 
quired, it  is  assumed  thut  the  modification  is  designed  to 
the  element’s  original  predicted  average  fatigue  life.  This 
function  indicates  whether  a modification  is  fatigue  tested 
before  its  incorporation  into  the  fleet.  It  is  an  input 
parameter  for  each  element  type, 

2 • 4 Actual  Average  Modified  Fatigue  life 

This  function  repeats  the  logic  established  in  Block 
1.5.  It  is  based  on  the  assumption  that  the  design  analysis 
of  the  modification  is  similar  to  that  of  the  original 
element  but  has  an  increased  probability  of  being  accurate. 
Therefore,  the  analysis  for  the  modification  results  in  an 
actual  average  fatigue  life  closer  to  the  required  life 
because  of  experience  gained  during  operational  usage.  This 
increased  probability  of  being  accurate  is  accounted  for  in 
SAIF  Is  by  decreasing  the  standard  deviation  and  increasing 
the  mean  of  the  log-normally  distributed  correction  factor 
discussed  In  Section  1.5.  The  standard  deviation  is  de- 
creased by  151  and  the  mean  is  Increased  by  15“  of  the  quan- 
tity (1.0  - mean), 

Monte  Carlo  techniques  are  again  used  along  with  the 
distribution  established  In  Figure  5 to  determine  the  actual 
average  fatigue  life  of  the  modified  element. 

2 • 5 Actual  Average  Modified  Fatigue  Life  Set  Hqual  to 
Free  Ic  t cTTTTTc 

If  the  result  of  Block  2.5  was  to  fatigue  tost  the 
modified  element,  then  it  Is  assumed  that  the  modified 
element  will  attain  its  predicted  life  or  be  redesigned  and 
retested  until  it  does.  Therefore,  the  average  lift*  pre- 
dicted by  Block  1.4  becomes  the  actuul  average  fatigue  lifo 
of  the  modified  element.  The  logic  then  returns  to  Block 
1.8  where  a fatigue  life  is  assigned  to  each  modified  ele- 
ment as  the  modification  is  installed  on  each  aircraft.  As 
was  done  for  the  original  clement,  the  individual  fatigue 
lives  are  determined  by  using  a Monte  Carlo  technique. 
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3. 


Derailed  Description  of  Block  3.0,  Reduce  Fatigue  hi  fe 
Becausc""6T  Production,  Sc rvj e e , cny  (lorrosH mi  Defects 
(,scc~  Figure  1 (ij  ' 


Figure  10.  Detailed  Flow  Diagram  - Reduce  Fatigue  Life 

Because  of  Production,  Service,  or  Corrosion  Defect 


3.1  Element  Has  Production  Defect 


As  one-time  occurrences,  production  defects  result  in 
structural  damage  only  when  thoy  initiate  the  progressive 
fatigue  failure  mechanism,  Typical  production  errors  in- 
clude surface  irregularities,  such  as  burrs,  nicks,  and 
gouges;  incorrect  dimensions  and  dimensional  tolerances; 
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improper  surface  finish  and  heat  treatment;  and  missing  or 
improperly  installed  fasteners  and  shims.  The  probability 
of  a given  element  having  a production  defect  before  enter- 
ing service  is  assumed  to  he  const  an  1 for  all  elements  of  a 
particular  structural  type  (spars,  frames,  etc,)  regardless 
of  individual  aircraft  or  the  aircraft  type. 

Data  establishing  the  production  defect  probability  for 
each  structural  element  type  were  obtained  from  the  Mechan- 
ical Reliability  Reports  (MRU's)  and  Service  Difficulty 
Reports  (SDK's)  in  conjunction  with  the  service  bulletins 
issued  periodically  by  the  airframe  manufacturers  and  Air- 
worthiness Directives  issued  by  the  FAA.  Because  of  the 
limited  regulatory  requirements  imposed  on  the  air  carriers 
in  documenting  the  nature  of  a structural  defect,  the  de- 
scription of  a defect  occurrence  in  MRR/SDR  data  is  gen- 
erally less  thorough  than  the  description  of  a potential 
structural  failure  in  the  pertinent  service  bulletin. 

A reported  fatigue  crack  is  classified  as  having  been  in- 
duced by  a production  defect  only  if  a recognizable  pro- 
duction error  either  is  explicitly  stated  in  the  MRR/SDR 
data  to  he  the  cause  of  the  failure  or  is  implied  to  be  the 
cause  by  reference  to  the  applicable  service  bulletin  or 
Airworthiness  Directive. 

From  a review  of  the  MRR/SDR' s submitted  to  the  FAA 
during  the  period  1963  to  1973,  a total  of  59  fatigue  fail- 
ures were  identified  which  could  be  attributed  to  production 
defects  with  almost  all  of  those  failures  resulting  from 
missing  or  improperly  installed  fasteners  and  shims.  The 
number  of  production  defects  found  in  each  type  of  struc- 
tural element  throughout:  the  fleet,  along  with  the  average 
number  of  individual  elements  in  each  aircraft  type,  are 
presented  in  Table  2.  The  MRR/SDR's,  and  consequently 
Table  2,  consider  only  production  defects  that  result  in 
cracks  or  corrosion.  Those  defects  that  do  not  affect  the 
aircraft  structural  integrity  are  not  considered  to  be 
significant  in  the  SAIb'li  logic. 

The  probability  that  an  individual  structural  element 
enters  service  with  a production  defect  is  equal  to  the 
ratio  of  the  number  of  production  defect  occurrences  iden- 
tified in  that  element  type  to  the  total  number  of  elements 
of  that  typo  within  the  fleet.  The  equation  for  determining 
the  probability  of  a production  defect  occurring  is 

1 p " fjiviT ■ * no. ■ ' u -nniTrnrnraar iiromont « ( 1 


For  the  pmposos  of  tin  sc  calculations,  it  is  assumed  that 
the  fleet  size  is  140b  aircraft,  the  largest  number  of 
pertinent  aircraft  registered  to  certified  route  air  car- 
riers in  any  given  year  during  the  period  1963  to  1973. 


These  data  are  from  Reference  6,  the  PAA  Statistical  Hand- 
book, 1973.  As  shown  in  Table  2,  the  largest  number  of 
production  defects  were  found  in  the  fuselage  stringers, 
although  the  fuselage  main  frames  have  the  greatest  proba- 
bility of  entering  service  with  such  a defect.  Since  it  is 
assumed  that  any  given  element  has  a small  finite  probabil- 
ity of  a production  defect  occurrence,  the  smallest  calcu- 
lated probability  (1.19xlO"S)  is  assigned  to  those  struc- 
tural elements  in  which  no  production  defects  were  identi- 
fied during  the  review  of  MRR/SDR  data. 


TABLE  2.  PROBABILITY  OF  A PRODUCTION  DEFECT  OCCURRING 


No , of 

Ave.  Mo. 

Probability  of 

Production 

Elements 

Production  Defect 

Element  Tyne 

Defects 

Per  A/C 

in  Indiv.  Element 

Fuselage 

Door  frame 

0 

10 

1.19  x 1 0 ' 1 * 

Window  frame 

0 

50 

1.19  x 10-'* 

Main  frame 

17 

180 

6.7  2 X 10- ! 

Floor  beam 

• 1 

60 

1.19  x 1 0 • s 

Keel  beam 

2 

60 

2,37  x 10‘s 

Pressure  web 

4 

6 0 

4,74  x It)’* 

Stringer 

24 

180 

9,48  x 10" 5 

Wing 

Access  frame 

0 

50 

1. 19  x 10"!* 

Rib 

2 

100 

1.42  x 10"! 

Spar 

3 

100 

2,13  x 1 0 " » 

Stringers 

3 

100 

2,13  x UK* 

Wing  Center  Section 

Rib 

1 

21 

3,39  x 10- # 

Spanwise  beam 

0 

21 

1,19  x 1 0 " s * 

Stringer 

1 

21 

3, 39  x 10" 6 

* estimated 


By  using  Monte  Carlo  techniques,  each  individual  struc 
tural  element  is  tested  at  the  time  an  aircraft  enters 
service  to  determine  whether  a production  defect  is  present 
If  tho  uniformly  distributed  random  number  drawn  is  less 
than  or  equal  to  the  appropriate  production  defect  occur- 
rence probability,  the  element  is  said  to  have  a production 
defect;  otherwise,  the  element  is  assumed  to  be  free  of  sue 
defects . 

3 . 2 Actual  Fatigue  Life  Reduced 

As  previously  stated,  a production  defect  results  in 
structural  damage  only  when  it  initiates  the  progressive 


fatigue  failure  mechanism,  Fatigue  crack  initiation  result- 
ing from  a production  defect  generally  occurs  early  in  the 
life  of  an  aircraft.  for  the  purposes  of  tins  simulation, 
it  is  assumed  that  the  effect  of  a production  defect  is  to 
lower  the  actual  fatigue  life  of  the  element. 

To  determine  the  fatigue  life  of  a production-damaged 
element,  the  distribution  of  fatigue  crack  occurrences 
resulting  from  production  defects  with  respect  to  time  of 
crack  initiation  was  determined.  With  the  use  of  two  as- 
sumptions regarding  crack  growth  rate  und  crack  size  at 
detection,  this  distribution  may  be  approximated  from  the 
available  MRR/SDR  data. 

Of  the  59  fatigue  cracks  attributed  to  production 
defects,  only  30  arc  documented  with  both  the  measured  crack 
length  and  the  aircraft  flight  hours.  Since  for  three 
fatigue  cracks  of  known  length  no  detection  time  was  re- 
ported, a time  was  assigned  to  each  on  the  basis  of  the 
aircraft's  year  of  manufacture,  the  date  of  the  crack  detec- 
tion, and  an  assumed  yearly  average  of  3000  flight  hours. 
Since  the  crack  length  was  not  reported  for  the  remaining  26 
fatigue  defects,  it  was  assumed  thut  the  distribution  of 
lengths  for  those  26  cracks  was  the  same  as  that  reported 
for  the  33  cracks. 

For  each  of  the  59  fatigue  cracks  initiated  by  produc- 
tion defects,  the  crack  length  is  plotted  as  a function  of 
the  crack  detection  time  in  figure  11.  Then  a crack  growth 
rate  was  postulated  by  constructing  a line  between  the 
origin  and  one  of  the  data  points  such  that  parallel  growth 
rate  lines  passing  through  each  of  the  other  data  points 
yield  no  negative  times  to  crack  initiation. 

Assuming  that  this  crack  growth  rute  is  constant  for 
all  fatigue  failures  initiated  by  production  defects,  the 
time  to  crack  initiation  can  he  computed  from  the  reported 
crack  detection  time  and  crack  length  by  the  following 
equation : 


1)1:.T 


+ (.log  0 . I - log 


£CU1  KCR 


(18) 


where 


KCI<  ■ 
KCI<  ■ 


time  to  crack  initiation  (flight  hours) 
crack  detection  time  (flight  hours) 
crack  length  at  detection  (inches) 
cruck  growth  constant 
1575  for  production  damaged  elements 
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figure  12  is  a histogram  of  times  to  crack  initiation 
determined  by  applying  liquation  (18)  to  each  of  the  Si) 
fatigue  cracks  initiated  by  production  defects.  Several 
standard  statistical  models  were  fitted  to  the  data  in 
figure  12,  with  the  Weibull  distribution  offering  the  best 
fit.  Thus,  whenever  the  simulation  program  determines  that 
a particular  element  has  a production  defect,  the  time  to 
first  crack  initiation  is  drawn  from  the  above  Weibull  dis- 
tribution rather  than  from  the  original  Weibull  distribution 
discussed  in  Section  1.8.  Times  to  second  and  third  crack 
initiations  are  still  drawn  from  the  Weibull  distribution  in 
Section  1.8. 


figure  17.  Histogram  of  Crack  Occurrences  on  Production 
Damaged  Elements 


3.3  iileincnt  Incurs  Service  Damage,  fatigue  Life  Reduced 

Like  production  defects,  service  defects  result  in 
structural  degradation  because  they  initiate  the  progressive 
fatigue  failure  mechanism.  Typical  examples  of  service 
damage  include  defects  occurring  during  normal  ground 
service  operations,  such  as  tears  and  dents  in  the  air- 
craft's skin  or  cargo  floor,  and  defects  occurring  durinb 
regular  maintenance  operations,  such  as  damage  to  parts 
during  installation  or  removal.  Data  available  from  MRR/ 
SDK's  show  that  the  service  damage  occurrence  rate  is 
constant  over  the  life  of  the  aircraft. 

A fatigue  crack  is  classified  as  having  been  induced  by 
service  damage  only  if  the  MRR/SDR  report  states  the  cause 


I to  be  a recognizable  service  damage  defect.  Of  the  61  crack 

defects  identified  that  could  be  attributed  to  service 
damage,  only  23  are  documented  with  both  the  measured  crack 
length  and  the  aircraft  flight  hours.  No  crack  detection 

* time  was  reported  for  five  fatigue  cracks  of  known  length, 

; and  the  measured  crack  length  was  not  reported  for  32  other 

f fatigue  defects;  for  the  one  remaining  occurrence,  neither 

i;  the  crack  length  nor  the  crack  detection  time  was  recorded. 

!■'  As  was  done  for  production  defects,  a crack  detection  time, 

l based  on  the  aircraft's  time  in  service  and  an  assumed 

[ average  number  of  flight  hours  per  year,  was  assigned  as 

l required.  In  addition,  it  was  assumed  that  the  distribution 

of  lengths  for  the  32  reports  where  crack  length  was  not 
;■  recorded  was  the  same  as  that  for  the  23  cracks  where  the 

| length  was  reported. 

I 

For  each  of  the  cracks  initiated  by  service  damage,  the 
I;  crack  length  is  plotted  in  Figure  .13  ns  a function  of  the 

crack  detection  time.  A crack  growth  rate  is  then  postu- 
lated such  that  no  negative  times  to  crack  Initiation  re- 
sult. Assuming  that  this  crack  growth  rate  is  the  same  for 
all  fatigue  failures  initiated  by  service  damage,  the  time 
to  crack  initiation  can  he  computed  from  the  reported  crack 

* detection  time  and  crack  length  by  the  following  equation: 
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where  ly^  a 1430  for  service-damaged  elements 

I: 

l- 

Assuming  that  service  damage  and  crack  initiation 
occur  almost  simultaneously,  the  times  determined  from 
liquation  (19)  can  be  used  to  construct  a histogram  of 
times  to  service  damage.  As  indicated  in  Figure  14,  the 
service  damage  occurrence  rate  is  independent  of  air- 
craft service  time. 

The  probability  of  service  damage  occurring  on  an 
element  of  a given  type  is  equal  to  the  number  of  occur- 
rences recorded  in  MRR/SDR's  divided  by  the  product  of  the 
number  of  flight  hours  in  the  data  and  the  number  of  indi- 
vidual elements  in  each  aircraft.  The  equation  for  calcu- 
lating the  service  damage  probability  is 

No.  of  Defects  * ly 

= 43,791,111  * No.  of  Individual  F. lemon ts 


where  ly  - adjustment  factor  discussed  below 

45,791,114  - total  flight  hours  in  MRR/SRR  data  base 

2 9 
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Figure  14.  Histogram  of  Service  Damage  Occurrences 


The  appropriate  occurrence  rate  for  each  structural 
element  type  along  with  the  number  of  occurrences  and  the 
number  of  individual  elements  In  each  aircraft  Is  given  in 
Table  3.  The  service  damage  occurrence  rate  is  adjusted  by 
a factor  of  two  because  of  the  results  of  the  Maintenance 
inspectors  Survey  (Volume  111)  which  indicates  that  service 
damage  is  twice  as  prevalent  as  actually  reported.  The 
difference  between  the  survey  results  and  service  damage 
reported  by  MRR/SDR  is  attributed  to  the  .limited  defect 
description  in  the  MRU/. SDR  reporting  format.  Therefore',  the 
MRR/SDR  data  are  used  to  determine  the  relative  rate  of  ser- 
vice damage)  occurrence  between  the  basic  element  types,  and 
the  survey  data  are  used  to  adjust  these  rates  to  what  Is 
felt,  to  he  more  realistic  occurrence  rates. 


No  fatigue  cracks  attributable  to  service  damage  are  1 
reported  for  four  element  types.  It  is  assumed  that  the  I 
wing  center  section  ribs  and  stringers  do  not  experience  j 
service  damage.  Although  no  service  damage  occurrences  in  I 
the  fuselage  pressure  webs  or  wing  access  frames  were  docu-  I 
mented,  it  is  assumed  that  individual  elements  of  these  1 
types  have  a finite  probability  of  experiencing  service  J 
damage  during  their  lives  and  that  the  occurrence  rates  for  I 
the  foregoing  are  the  same  as  those  for  the  fuselage  window  J 
frames  and  wing  stringers,  respectively.  1 

Using  Monte  tiarlo  techniques,  each  individual  struc-  | 
tural  element,  is  tested  when  an  aircraft  enters  service  to  J 
determine  whether  service  damage  is  incurred  and,  if  so,  at  I 
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what  time  the  resultant  fatigue  crack  initiates.  The  gov- 
erning equation  which,  in  conjunction  with  the  selected 
uniformly  distributed  random  number,  determines  this  time  to 
crack  initiation  is 


t - i ln(RN) 

A 


(21) 


where  X is  the  appropriate  service  damage  occurrence  rate 
for  the  element  type,  and  RN  is  the  random  number  selected 
by  Monte  Carlo  methods. 


TABLE  3.  SERVICE  DAMAGE  OCCURRENCE  RATES 


Element  Type 

No.  of  Service 
Damage  Occurrences 
Before  11040  Pit  Hr 

Avo.  No. 
moments 
Per  A/C 

Ad.j  . Service 

Damage  Occurrence 
Rate  focc./flt.brl 

Puseliiye 

Door  frame 

7 

10 

3.06 

X 

10'  * 

Window  frame 

8 

SO 

6.09 

X 

in-  • 

Main  frame 

21  . 

180 

5.10 

X 

io- 1 

floor  beam 

1 

60 

2.4  3 

X 

10  1 “ 

Keel  beam 

1 

60 

2.43 

X 

10  ■ 1 0 

Pressure  web 

0 

60 

2.4  3 

X 

10'  1 5* 

Str.inyor 

14 

180 

3.4  0 

X 

10-  1 

Wing 

Access  frame 

0 

50 

1.31 

X 

10-’* 

Rib 

1 

100 

4.36 

X 

10-  1 0 

Spar 

1 

100 

4.36 

X 

in- 1 c 

Stringer 

3 

100 

1.31 

X 

10'" 

Winy  Center  Section 

Rib 

0 

21 

0 

Spanwiso  beam 

4 

21 

M 

to 

CD 

X 

10*  * 

Stringer 

0 

21 

0 

* estimated 


3 . 4 Corrosion  Initiation  Occurs 

Each  individual  structural  element  has  a finite  prob- 
ability of  experiencing  corrosion  during  its  service  life. 
The  corrosion  occurrence  rate  for  each  element  type  is 
determined  from  the  service  experience  documented  in  the 
MRR/SDR's.  The  first  step  in  the  formulation  of  the  corro- 
sion occurrence  rate  is  a determination  of  the  corrosion 
growth  rate.  For  this  determination,  the  corrosion  occur- 
rence data  of  the  DC-9  wing  center  section  stringers  are 
analyzed  by  examining  the  reported  corrosion  depths  as  a 
function  of  the  aircraft's  accrued  flight  hours.  Of  the  52 
reported  corrosion  occurrences  on  the  DC-9,  only  36  are 


' 
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\ 
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documented  witli  both  the  corrosion  depth  and  aircraft  life. 
Furthermore,  in  15  of  these  occurrences,  the  corrosion 
initiated  at  fastener  holes  in  the  stringers.  Since  the 
corrosion  depth  growth  rate  may  he  related  to  the  state  of 
stress,  the  instances  in  which  corrosion  occurred  in  a 
stress  concentration  are  excluded  during  the  initial  analy- 
sis. 


For  each  of  the  remaining  .11  corrosion  occurrences,  the 
reported  corrosion  depth  is  plotted  as  a function  of  the 
corrosion  detection  time  (sec  Figure  I 5 j . A corros  ion 
growth  rate  is  then  postulated  hy  const  met i ng  a line  be- 
tween the  origin  and  one  of  the  data  points  such  that  paral- 
lel growth  rate  lines  passing  through  each  of  the  other  data 
points  yield  no  negative  times  to  corrosion  initiation. 

This  is  a somewhat  conservative  approach  in  that  it  assumes 
corrosion  may  initiate  as  soon  as  an  aircraft  enters  ser- 
vice. The  calculated  rate  is  1.0!)  x 10*5  inches  per  flight- 
hour  . 


Figure  15.  Time  of  Defect  Detection  vs.  Depth  of  Corrosion 
for  Corrosion-Damaged  lilements 
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Now  the  time  to  corrosion  initiation  may  be  determined 
each  occurrence  of  known  depth  and  detection  time  From 
following  equation: 


= t 


co 


CO 


CO 


(22  j 


whore  t0  » time  to  corrosion  initiation  (flight  hours) 
t = corrosion  detection  time  (flight  hours) 
d,0  « corrosion  depth  at  detection  (Inches) 
k ■ corrosion  growth  rate  (inches  per  flight  hour) 


for  those  corrosion  occurrences  of  unreported  detection 
time,  a value  of  tco  is  assigned  on  the  basis  of  the  air- 
craft's time  in  service  and  an  assumed  average  number  of 
flight  hours  per  year,  l!or  those  corrosion  occurrences  of 
unrecorded  depth,  a value  equal  to  the  mean  of  the  21  pre- 
viously plotted  corrosion  occurrences  is  assigned.  The 
corrosion  Initiation  time  for  each  of  the  other  reported 
corrosion  occurrences  throughout  the  fleet  can  now  be  cal- 
culated by  using  the  same  equation  and  growth  rate. 

It  should  be  noted  that  the  same  corrosion  depth  growth 
rate  is  used  to  determine  the  time  to  corrosion  initiation 
for  all  occurrences  regardless  of  whether  the  corrosion  is 
located  in  a stress  concentration.  This  procedure  was  found 
acceptable  after  analysing  the  data  for  corrosion  coincident 
with  stress  concentrations , Tor  the  IS  fully  documented 
occurrences  Ln  the  HC-9  wing  center  section  stringers,  t ho 
previously  established  techniques  wore  then  applied  to  those 
data  to  postulate  a new  corrosion  depth  growth  rate.  This 
rate  was  found  to  be  very  close  to  1 . 09  x 10*5  in.  per 
flight  hour.  Since  this  rate  was  similar  to  that  previously 
calculated,  the  original  growth  rate  was  chosen  to  conserva- 
tively define  time  to  corrosion  initiation  for  all  occur- 
rences . 


l-'or  each  structural  element  type,  a curve  of  cumulative 
corrosion  occurrences  throughout  the  fleet  is  plotted  as  a 
function  of  time  to  corrosion  initiation,  linen  plot  ex- 
hibits cither  a constant  corrosion  occurrence  rate  or  two 
constant  but  unequal  rates.  The  occurrence  rate  during  the 
initial  service  life  of  the  aircraft  is  frequently  lower  be- 
cause protective  finishes  and  coatings  are  Intact  and  have 
not  deteriorated  or  been  scratched  and  gouged.  As  the 
protective  finishes  break  down,  the  rate  increases.  At  the 
nigh-occurrence  rates,  the  data  shows  that  a time  is  reached 
when  the  number  of  high- time  aircraft  in  the  fleet  decreases 


and  accordingly  the  occurrence  rate  also  decreases.  This 
decrease  is  ignored  in  the  simulation  since  all  the  aircraft 
in  the  fleet  being  simulated  are  flown  to  their  retirement 
lives.  These  corrosion  occurrence  curves  are  shown  in 
Figures  16  through  27.  The  occurrence  rate  for  a fleet,  of 
1406  aircraft,  is  then  derived  from  this  data  by  using  liquation 
(23).  The  fleet  size  data  were  obtained  from  the  h'AA  statis- 
tical handbook  of  aviation  for  1070  and  1973  (References  6 
and  7). 

„ A Occurrences  . , 
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3,5  Corrosion  Occurs  in  Stress  Concentration 

The  fatigue  life  reduction  of  a structural  element 
resulting  from  corrosion  damage  depends  upon  the  stare  of 
stress  in  tho  cn'-ruded  region,  for  all  elements  of  a 
particular  structural,  type,  the  probability  of  corrosion 
occurring  within  a stress  concentration  is  assumed  to  be 
constant  over  the  life  of  the  aircraft,  regardless  of  indi- 
vidual aircraft  or  aircraft  typo. 
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Figure  16.  Corrosion  Occurrence  Rate  for  fuselage  floor  Hearns 
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Figure  21.  Corrosion  Occurrence  Rate  for  luse  I ago  Stringer 
(Hot  tom) 
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Corrosion  Occurrence  Rate  for  Fuse  I a ge  St  r in  iter 
(Side  and  Top) 


Figure  22 
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Figure  27.  Corrosion  Occurrence  Rale  for  Wing  Stringers 


The  corrosion  occurrences  documen  t cd  in  i lie  MKK/SDk's 
are  class  If  Led  as  being  located  within  a stress  concentra- 
tion only  If  a recognisable  design  feature  known  to  he  a 
stress  riser  is  coincident  with  the  reported  corrosion.  The 
two  types  of  stress  riser  identified  in  the  available  data 
"o  fastener  holes  and  bend  radii  lor  fillets). 

Of  the  583  documented  corrosion  occurrences,  75  were 
classified  its  being  located  in  stress  concentrations.  The 
number  of  corrosion  occurrences  within  stress  concentrations 
found  in  each  type  of  structural  element  as  well  as  the 
total  number  of  corrosion  occurrences  in  elements  of  (hat- 
type  are  presented  in  Table  I.  The  probability  that  corro- 
sion exists  in  a stress  concentration  is  equal  to  the  ratio 
of  the  number  of  corrosion  occurrences  found  in  stress 
concentrations  to  the  total  number  of  corrosion  occurrences 
identified  in  that,  structural  element  type  within  the 
fleet.  Since  it  is  assumed  that  a corrosion  occurrence  in 
any  typo  of  structural  element  has  a finite  probability  of 
being  located  in  a stress  concentration,  a probability  must 
l)e  assigned  to  those  structural  elements  for  which  no  oc- 
currences of  corrosion  are  reported.  This  assignation  is 
made  by  assuming  a probability  equivalence  between  element- 
types  of  generally  similar  design  and  structural  function. 
For  example,  the  probability  of  a corrosion  occurrence  being 
located  in  a stress  concent  rat i on  of  either  a door  frame  or 


a window  L'rame  is  assumed  to  be  t lie  same.  Probability 
values  assigned  in  this  manner  are  indicated  by  an  asterisk 
i.  n T a b .1.  e 4 . 


TABU:  4.  PROBABILITY  OF  CORROSION  OCCURRENCE  IN  A 
STRESS  CONCENTRATION 


Structura 1 
Element  Typo 

No.  of 
Corrosion 
Occurrences 

No . of 

Occurrences  in 
Stress  Concen. 

fuselage 

• 

Door  frames 

] 3 

5 

Main  framer. 

23 

3 

Window  frames 

3 

2 

Floor  beams 

53 

1 

Keel  beams 

5 

0 

Pressure  wehs 

1 

0 

Stringers 

257 

35 

W ing 

Access  frames 

0 

0 

Itibs 

9 

1 

Spars 

57 

3 

Stringers 

46 

3 

Wing  Center  Section 

Ribs 

0 

0 

Spanwisw  beams 

18 

1 

Stringers 

1,2 

2 1 

Probability  of 
Corrosion  Occur, 
in  Stress  Concen, 


0 . 3 fl  5 
0.130 
0 . 3H5* 
0 . U 1 7 
0.U17* 
0. 136* 
0 . 136 


0.385* 

n.ll  I 
0.053 
0 . 06  5 


0 .036* 

0.056 

0.238 


* estimated 


Using  Monte  Carlo  methods,  each  incident  of  corrosion 
Is  tested  at  the  time  it  occurs  to  determine  whether  it  is 
located  in  a stress  concentration.  if  the  uniformly  dis- 
tributed random  number  drawn  is  less  than  or  equal  to  the 
appropriate  probability  of  corrosion  occurrence  in  a stress 
concentration,  the  corrosion  is  assumed  to  be  located  in  a 
stress  concentration;  otherwise,  the  corrosion  is  assumed  to 
occur  in  a uniform  stress  field. 

3 . 6 Actual  Fatigue  Life  Reduced  Because  of  Corrosion  in  a 

5 1 r e s s'  C o n c e n~t~r~a  t i o n 

The  presence  of  corrosion  on  a structural  element 
contributes  to  the  potential  failure  of  the  element  by 
reducing  the  original  fatigue  life  of  the  element.  Tests 
conducted  on  spar  caps  taken  from  HU- 16  aircraft  and  docu- 
mented in  Reference  8 show  that  severe  exfoliation  corrosion 
reduces  fatigue  life  significantly,  but  that  surface  pitting 
or  very  mild  exfoliation  has  only  a minor  effect  on  fatigue 
life. 
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When  it  has  born  determined  that  cur  rosiiin  has  nee  u r ml 
in  a stress  concentration,  the  fatigue  life  of  the  element 
is  reduced  by  a factor  that  is  an  input  variable  in  Block 
1.1.  Reference  ts  shows  that  when  corrosion  occurs  in  a 
stress  concentration,  it  generally  manifests  itself  as 
exfoliation.  The  reference  indicates  that  fatigue  tests 
conducted  on  MU-16  spar  caps  showed  that  severe  exfoliation 
corrosion  results  in  a fatigue  l ife  reduction  of  up  to  7 . 

However,  since  the  tests  wore  conducted  on  specimens  that 
had  previously  experienced  service  fatigue  damage,  it  was 
felt  that  approximately  3()'i  of  the  reduction  was  because  of 
such  damage.  Therefore,  40'i,  of  the  fatigue  life  reduction 
was  attributed  to  exfoliation  corrosion  and  the  suggested 
simulation  input  for  fatigue  life  reduction  when  corrosion 
occurs  in  a stress  concentration  is  0.40. 

3 . V Actual  l; a t ign c hi  f e _ Reduced  lie  c a u s e of  Co r ro s ion 

Outside  a Stress  iTohcenT  ratio  if 

Corrosion  outside  a stress  concentration  also  a f fee  ts 
the  fatigue  life  of  the  element,  hut  the  fatigue  life  reduc- 
tion is  less  severe  than  when  the  corrosion  is  in  u stress 
concentration.  Reference  1 indicates  that  the  fatigue  life 
reduction  because  of  corrosion  outside  a stress  concentra- 
tion is  approximately  one  I If  of  the  fatigue  life  reduction 
when  the  corrosion  is  in  a stress  concent  rat i on . Therefore, 
the  suggested  simulation  input  for  fatigue  life  reduction 
when  corrosion  is  outside  a stress  concentration  is  0 . It) . 


r*  rao i RiJW-TKC^crwtTJir ' Km.^^SFcrmRfTfTSfii  rcs^raraff 


4 . Detailed  Desc  r Lpt  ion  of  Block  4_.  0 .Reduce  Strength  Because 
of  Crack'  Crowth/T’f  ojee  t T fie  t o Pa  1 l ure  (see  1;T  g u f e 2 8 ) 


Figure  28.  Detailed  Flow  Diagram  - Reduce  Strength  Because 

of  Crack/Corrosion  firowth/Predict  Time  to  Failure 


4 . 1 Predicted  Corrosion  Growth  Rate 

If,  as  a result  of  Block  .".4,  it  has  been  determined 
that  corrosion  is  present  on  an  element,  this  function 


44 


~*il 


defines  the  rate  at  which  the  area  and  depth  of  the  cor- 
rosion will  grow.  The  growth  rates  presented  are  the  result 
of  an  analysis  of  MIIR/’S1)U  data.  Because  of  the  small  number 
of  usable  dutu  points,  the  growth  rates  presented  are  de- 
termined from  data  for  different  element  typos  and  are 
assumed  to  represent  a nominal  condition.  This  nominal 
growth  rate  is  modified  by  applying  u correction  factor  to 
each  element  type  based  on  the  corrosion  resistance  rating 
(CRR)  as  determined  by  the  MSG- 2 conference  convened  for  the 
certification  of  each  aircraft  type.  The  CRR  has  a range  of 
1 through  4 with  4 representing  the  elements  that  are  most 
resistant  to  corrosion.  The  growth  rate  for  each  element 
then  becomes 

Actual  Rate  - Nominal  Rate  x Adjustment  factor  (24) 

where 

CRR  Adjustment  Fac  t or 

J 1.50 

2 1.2  5 

3 1.00 

4 0.75 

Nominal  depth  rate  B 1.09  x 10*5  In.  /.fit  hr 
Nominal,  area  rate  “ 2.0  x .10*3  sq.ln./flt  hr 


The  corrosion  depth  growth  rate  was  derived  in  Block 
3.4  from  the  data  presente.,  in  Figure  15.  The  corrosion 
area  growth  rate  was  also  determined  from  the  MRR/SDR  data. 

4 . 2 Average  Crack  Growth  Rato 

The  simulation  program  uses  two  constant  crack  growth 
rates  to  approximate  the  normally  non-constant  fatigue  crack 
growth  rate.  Those  two  constant  rates  are  inputs  for  each 
element  type.  The  first  constant  rate  represents  the  crack 
growth  from  crack  initiation  until  the  critical  crack  length 
is  reached.  This  is  the  slow  growth  rate  period. 


The  second  constant  rate  represents  the  fast  growth 
rate  and  accounts  for  the  crack  growth  from  critical  crack 
length  to  structural  failure.  Taken  from  References  9 and  10, 
Figures  29  through  3a  are  typical  crack  growth  curves  for 
selected  components  of  a DC-8  and  DC-10.  Drawn  over  the 
crack  growth  curves  in  each  of  the  figures  are  the  two 
constant  approximations  used  in  the  SAl'Fli  program.  For 
those  very  few  elements  (Figure  35)  whose  crack  growth 
curves  do  not  conform  to  the  general  shape  discussed  above, 
a single  constant  approximation  is  used.  To  determine  the 
input  parameters  required  for  the  simulation  program,  the 
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slow  growth  portion  of  the  crack  growth  curve  must  he  ex- 
trapolated backward  to  zero  crack  length  to  determine  the 
constant  rate  from  crack  initiation  to  critical  crack 
length.  The  fast  growth  portion  of  the  crack  growth  curve 
starts  at  the  critical  crack  length  and  is  applied  until 
structural  failure. 
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FAIL  SAIL 

PANEL  DEVELOPMENT  TESTS  AND 
PROJECTED  AIRCRAFT  UTILIZATION 
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Figure  29.  Two-Bay  Longitudinal  Crack  Propagation  Rate  for 
Typical  Lower  Surface  Wing  Stringers 
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Figure  30.  Crack  Propagation  Rule  for  Typical  Fuselag 
Main  Frames 
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Figure  31.  Crack  Propagation  Rate  for  Typical  Fuselage 
Stringers 
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Figure  32.  Crack  Propagation  Rate  Cor  Typical  Lower 
Rear  Wing  Spars 


Figure 


4.3  Actual  liloinent  Crack  Growth  Rate 


Crack  growth  rate,  like  fatigue  life,  is  a probabilis- 
tic variable  and  must  be  treated  statistically.  The  scatter 
present  in  the  growth  rate  distribution  reflects  both  the 
Dftsic  characteristics  of  the  fatigue  process,  the  corrosion 
environment  variation,  and  the  load  environment  variation  of 
typical  aircraft  structures.  On  the  basis  of  the  results 
from  full-scale  testing,  liggwertz  (Reference  11)  has  deter- 
mined that  the  standard  deviation  for  the  crack  growth  rate 
is  approximately  one-half  that  of  the  fatigue  life.  Taken 
from  Reference  12,  Figure  37  characterizes  the  variability 
of  crack  initiation  and  growth  in  distributional  forms.  In 
the  simulation  both  the  slow  and  fast  crack  growth  rates  are 
assumed  to  be  normally  distributed  with  moans  equal  to  the 
slow  and  fast  growth  rates  determined  in  Section  4.2.  The 
standard  deviation  is  set  equal  to  one  half  that  used  in  the 
time  to  crack  initiation  distribution.  Both  of  these  param- 
eters. the  mean  slow  and  mean  Past  growth  rates,  are  re- 
quired as  SAlFli  input. 
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Figure  37.  Typical  Crack  Initiation  and  Growth  Rate 
Distribut ions 


f 


b; 

I 


if- 

I 


The  slow  and  fast  crack  growth  rates  for  each  indi- 
vidual element  on  each  aircraft  in  the  fleet  are  selected  by 
a Monte  Carlo  method.  The  selection  is  controlled  so  that 
the  fatigue  life,  slow  growth  rate,  and  fast  growth  rate  of 
an  individual  element  all  reflect,  a consistent  material  scat- 
ter and  load  environment.. 


4 , 4 Corrosion  Iff  foots  on  Strength 


References  8,  13,  and  14  present  test  results  whLch 
show  that  the  reduction  In  static  strength  because  of  cor- 
rosion is  negligible  until  the  loss  of  cross- sect Lona 1 area 
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becomes  an  extremely  significant  portion  of  the  element- 
cross  section.  A paper  present  oil  at  the  1072  T r i He  r v i ce 
Conference  on  Corrosion  (Reference  1 5 I present  eil  test  data 
on  the  effect;  of  corrosion  on  static  strength,  'I  lie  data 
scatter  was  so  large  that  it  must  be  concluded  that  the 
effect  of  corrosion  on  static  strength  cannot  be  measured 
accurately  by  present  standards. 

An  examination  of  the  growth  rates  defined  in  Block  <)  . 1 
and  the  detection  probabilities  defined  in  Section  5.0  indi- 
cates that  corrosion  will  certainly  be  detected  before  any 
detectable  reduction  in  static  strength  lias  taken  place. 
Reference  9 supports  this  finding  with  regard  to  detection. 

Reference  lo  examines  the  effects  of  corrosive  environ- 
ments on  fatigue  life  of  aluminum  alloys  under  maneuver 
spectrum  loading.  The  test  results  showed  that  in  general 
all  the  aluminum  alloy  plate  materials  tested  experienced 
significant  and  progressive  reductions  of  mean  fatigue  lives 
for  increasingly  severe  corrosive  environments.  The  average 
crack  propagation  rate  was  approximately  tripled  by  the 
static  and  cyclic  corrosion  environments.  However,  the 
effective  crack  length  as  measured  after  specimen  failure 
appeared  to  he  unaffected  by  environment.  This  suggests 
that,  the  effects  on  residual  cracked  strength  of  the  various 
corrosive  environments  is  negligible.  Therefore,  although 
it  may  be  intuitively  felt  that  corrosion  must  have  an 
effect  on  static  strength,  the  data  presently  available  does 
not  support  that  opinion.  It  Is  also  apparent  that  state- 
of-the-art  material  selection  and  preventive  coating  appli- 
cations have  eliminated  corrosion  as  a major  factor  In 
' ,i  tustroph  i c ace  i dents . 

While  corrosion  dons  noi  affect  the  relationship  be- 
tween crack  length  and  residual  strength,  it  does  accelerate 
strength  degradation  by  its  effect  on  crack  propagation 
rate.  When  corrosion  is  present,  the  crack  propagation  rate 
is  increased  by  the  same  factors  used  to  decrease  fatigue 
life  ns  defined  in  Sections  3.0  and  3.7. 


4.5  St  rung  tli  Reduction  Bocjuis_e_  uj’  Crack  ilrowth 

for  the  simulated  structural  elements,  it  is  assumed 
that,  the  original  ultimate  strength,  Su , is  constant  until 
the  time  of  crack  initiation,  t ° 0,  and  that  after  crack 
initiation  the  subsequent  residual  strength,  S,  ran  be 
expressed  as  a function  of  time  (flight  hours  1.  As  dis- 
cussed in  Section  1.2,  the  two-part  approximation  to  the 
growth  rate  of  a single  fatigue  crack  has  the  general  con- 
figuration shown  in  figure  38.  It  will  also  he  assumed  that 
the  relationship  between  crack  length  and  residual  strength, 
S,  is  as  shown  in  figure  39,  with  S going  to  1.0  at  the 
crack  length  corresponding  to  level  flight  failure.  If 
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the  curves  in  Figure  38  and  Figure  39  ure  combined , the 
resultant  composite  curve  illustrates  the  functional  rela- 
tionship between  time  (flight  hours)  and  S(t),  where  S(t) 
is  the  residual  strength  of  the  element  at  time  t ufter 
crack  initiation.  Such  a curve  for  a typical  element  is 
shown  in  Figure  40.  Approximately  10  percent  of  the 
structural  elements  under  consideration  have  no  change  in 
fatigue  crack  growth  rate  at  the  critical  crack  length. 

For  these  elements  the  functional  relationship  between 
flight  hours  and  residual  strength,  S(t),  is  illustrated  by 
the  two-part  curve  shown  in  Figure  41,  'Approximately  IS 
percent  of  the  structural  elements  under  consideration  have 
critical  crack  lengths  and  average  slow  crack  growth  rates 
such  that  neither  the  fail-safe  length  nor  the  critical 
crack  length  can  be  reached  in  the  aircraft's  service  life, 
For  these  elements  the  functional  relationship  between 
flight  hours  and  residual  strength,  S(t),  is  the  linear 
curve  shown  in  Figure  42. 
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If  there  is  a second  crack  initiation,  a two-part 
approximation  to  the  growth  rate  of  the  sum  of  the  two  crack 
lengths  is  calculated,  The  slow  growth  port  ion  of  the 
summed  crack  approximation  is  defined  hy  two  points:  The 
first  point  is  the  length  of  the  first  crack  when  the  second 
crack  initiates,  and  the  second  point  is  the  sum  of  the 
crack  lengths  when  the  first  crack  reaches  critical  crack 
length.  The  fast  growth  portion  of  the  approximation  is 
determined  hy  a Least  - .squares  fit  of  the  sum  of  the  crack 
lengths  when  the  second  crack  initiates,  when  the  second 
cruck  reaches  critical  length,  and  when  the  aircraft  is 
retired.  A new  S(t)  curve  is  now  determined  as  before, 
except  that  t*0  now  represents  the  time  of  the  second  crack 
initiation.  S('t»0)  is  now  equal  to  Su  minus  the  loss  of 
residual  strength  because  of  the  first,  crack  growth  before 
the  second  crack  initiation.  If  a third  crack  initiates, 
the  same  procedure  is  repeated. 

4 , 6 Load  llxcoeds  R e d u c o d S t_r e lyg t_h 

This  function  involves  a probabilistic  determination  uf 
the  maximum  flight  load  experienced  hy  an  airplane  and  the 
comparison  of  this  load  with  the  strength  of  the  elements  to 
project  time  to  failure  after  a crack  initiation.  The  data 
presented  in  Reference  1?  was  used  to  define  the  distribu- 
tion of  positive  and  negative  normal  accelerations  experi- 
enced hy  two  aircraft  equipped  with  NASA  VGII  recorders, 
which  provide  continuous  time-history  records  of  indicated 
airspeed,  normal  acceleration,  and  pressure  altitude,  The 
data  were  collected  over  a 2 -year  period  on  two  identical 
four-engine  turbojet:  transport  airplanes  during  routine 
commercial,  operations  of  a single  airline.  The  data  covered 
flights  mostly  over  the  eastern  half  of  the  Continental 
United  States  and  a few  to  the  West  Coast  and  to  northern 
South  America.  The  data  consisted  of  37(>(i  flight,  hours  ol 
operational  maneuver  and  gust  accelerations  anil  2 1 i) . 7 flight 
hours  of  check- flight  maneuver  accelerations.  These  data 
compared  closely  with  those  for  another  type  of  four  engine 
turbojet  transport,  The  operational,  maneuver,  operational 
gust,  and  check- f I ight  maneuver  accelerations,  both  positive 
and  negative,  were  combined,  and  the  exceedances  per  flight 
hour  were  calculated  for  each  deviation  from  level  flight  in 
O.lg  increments,  A 1 oust  - squn res  curve- fit  computer  program 
was  then  used  to  fit  an  exponential  curve  to  the  exceed- 
ances. The  equation  for  the  exponential,  curve  is  as  fol- 
lows: 


I’  ( S ) 15  A exp  | bS 

tl  il 


(25) 


where  1*  C Sf» ) is  the  number  of  flight  loads  per  hour  which 
exceed  the  load  level  S;1 , and  A and  b are  input  parameters. 
A plot  of  the  observed  exceedances  is  presented  in  figure 
43,  The  .least -squares  curve  fit  i adjusted  to  give  the 
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With  S (t)  as  expressed  above,  (lie  number  of  flight  loads  per 
hour  which  exceed  the  residual  strength  at.  time  t is 


P [>  ( t :]  “ A e x p 

’b(Su  ’ V1' 

* 

1 ■'  t , 

b [S(t)]  " A exp 

b S1  - iMt 

- V 

j 

,t  •-  t < t2 

P [s  ( t. )]  “ A e x p 

b Sf  - R - (.  v. 

• t2) 

,t  > t, 

Ld 

Adopting  the  same  assumptions  as  I.undbcrg  and  Hggwertz 
(Reference  18),  the  above  expressions  for  the  residual 
strength  exceedance  rate  can  be  substituted  for  the  risk 
function,  A ( t ) , in  the  reliability  formula.  Recall  fiom 
Volume  I tie  reliability  formula 


bit.')  ■=  1 - exp 


f A ( t. ) dt 


2 8 


Making  the  above  substitutions  vlelds 
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1 - exp 


c 2 m 


l’’(t)  “ 1 - oxp 


Jl  A exp  [ 


A exp  I bS  - bR^  t 


dt 


J A exp  - b R 2 C ^ “ tj)J 


dt 


(30) 


h K t 


F ( t ) - 1 * exp 


f2  A exp  j^bS^  - bR2(t  - t^jdt 


J 1 A oxp  £ bSu  - bll^tjdt 


(31) 


/ A exp  [bSf  • bRs(t  - 1 2 ) ]dt 


, t -■  t. 


where  P(t)  is  the  probubility  oi'  failure  before  time  t 

of  an  aircraft  structure  which  had  a crack  initia- 
tion at  time  t-0. 

To  calculate  the  predicted  times  to  failure  from  uni- 
formly distributed  random  numbers,  equate  the  random  number 
RN  to  F (t ) and  solve  for  t.  for  P(t)  as  expressed  above, 
and  setting  RN  equal  to  1 - F(t)  for  convenience,  solving 
for  t yields 


4 ■ ■ CT7  ln 


bRj  in  (RN) 
A exp  (b'Su) 


+ 1 


t «.  t. 
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bR. 


A exp  rbS1  + bR 2 1 


A exp  (bS,j 

,n(RN)  + ^ 


A exp  (bSu) 
BIT 


(exp  (-bR^tj.)  - l) 


, t1  < t < t.,  (33) 
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When  there  is  a second  crack  initiation,  a now  time  to 
failure  is  calculated  with  t *>  0 now  corresponding  to  the 
time  of  second  crack  initiation  and  Su  now  equaling  the 
ultimate  strength  minus  t he  loss  of  residual  strength  be- 
cause of  the  first  crack  growth.  The  same  procedure  is 
repeated  if  there  is  a third  crack  initiation. 
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Petal  led  Description  of  Mock  S . 0 , Periodic  Inspection 
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Figure  44.  Detailed  Flow  Diagram  - Periodic  Inspection  of 
Elements/Aircraft  Deleted  from  Fleet 
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5.1  Current  Periodic  Inspection  T nt_ct;va_l_s 

The  initial  inspection  intervals  are  input  parameters. 
The  inspection  intervals  can  be  those  recommended  by  the 
Maintenance  Review  Board  (MRB)  or  those  submitted  by  an  air 
carrier  as  part  of  the  Standard  Operations  Specification. 

The  simulation  is  designed  to  accept,  a standard  ATA  four- 
level  inspection  program.  The  four  levels,  A through  D,  arc 
defined  as  follows: 

A Preflight:  Visual  inspection  conducted  from 

ground  level  and  primarily  covering  lower  exterior 
wing  and  fuselage  surfaces. 

B Service:  Close  visual  inspection  of  bottom  of 
wing,  lower  fuselage,  top  of  wing,  and  known 
problem  areas.  This  also  may  include  the  front  of 
the  forward  spar  and  rear  of  the  aft  spar  in 
readily  accessible  areas. 

l!  Phase:  Close  visual  .inspection  of  aircraft  ex- 
terior and  easily  accessible  interior  areas,  such 
as  baggage  compartments  ami  door  frames.  NTH  of 
selected  areas  of  the  aircraft. 

D Overhaul:  Detailed  inspection  of  entire  aircraft. 
This  level  may  be  conducted  both  on  a sampling 
basis  and  during  several  separate  inspections. 


It  is  assumed  that  each  higher- level  inspection  Includes  all 
the  lower-level  inspections  down  to  the  lowest  level  speci- 
fied for  a particular  element  type.  Therefore,  if  a D-level 
inspection  is  being  conducted  and  the  lowest  interval  speci- 
fied for  that  element  was  a R- level,  then  the  current  In- 
spection would  include  the  B- level,  C-  level,  and  D-lev-  ! 
.inspections . 


Sampling  inspections  are  accounted  for  in  the  D- level 
inspection  logic  by  reducing  the  probability  of  defect 
detection.  The  amount,  of  reduction  is  directly  proportional 
to  tile  size  of  the  sample  inspected.  liach  time  an  aircrafr 
is  inspected  at  the  U- level,  the  probability  of  defect  de- 
tection is  multiplied  by  the  fractional  size  of  the  sample. 
This  decimal  fraction  is  an  input  parameter. 


The  J 
i n c 1 c a s e s 
depend i ng 
amount  of 
f req uency 
into rva  1 s 


ogle  that  follows  in  Block  7.1  automatically 
the  frequency  of  inspection  at  certain  levels 
upon  the  extent  of  the  defects  being  found.  The 
increase  is  an  input  function.  These  changes  in 
are  accounted  for  in  this  function,  and  the  now 
are  used  to  schedule  subsequent  inspections. 


b 1. 


5.2 


Next  Inspection  at  Overhaul  Level 

This  function  is  a check  of  the  inspection  level  being 
scheduled.  It.  is  present  because  the  installation  of  modi- 
fications is  normally  scheduled  only  during  overhaul  inspec- 
tions, and  the  modification  installation  takes  precedence 
over  the  inspection  process. 

5 . 3 Modification  J’cjnd  in  on  Thi  s JiJ  e m ojn 

Once  the  decision  is  made  to  develop  a modification,  it 
generally  takes  3 to  6 months  to  design  the  modification, 
procure  materials  and/or  parts,  and  set  up  the  tooling 
required  for  installation.  These  lead  times  are  input 
parameters  that,  depend  on  the  element  being  considered. 

Because  SAIFE  processes  time  in  flight  hours,  calendar 
days  or  months  must,  he  converted  to  the  equivalent  simula- 
tion time  by  liquation  (35)  or  (3b) 

Calendar  Hays  * 8.2  * equivalent  simulation  time  (35) 

Calendar  Months  * 250  * equivalent,  simulation  time  (36) 


If  the  modification  lead  time  requirements  are  satis- 
fied, the  logic  progresses  to  the  installation  routine,  and 
it  is  assumed  that  existing  defects  are  repaired  during  the 
modification  process.  If  the  lead  time  requirements  have 
not  been  satisfied,  or  if  no  decision  has  been  made  to 
modify  the  element,  the  logic  continues  through  the  inspec- 
tion subroutine. 


5 . 4 Defect  Internal  or  External 

The  Inspection  level  at  which  a defect  will  be  detected 
depends  partly  on  whether  the  defect  is  internal  or  exter- 
nal. Assuming  that  the  occurrence  of  erncks  or  corrosion  is 
a statistical  trial  and  that  each  trial  is  independent  of 
ail  other  trials,  then  the  probability  of  the  corrosion 
being  external  in  a single  occurrence  is  simply  the  number 
of  external  occurrences  divided  by  the  total  number  of 
occurrences.  Based  on  the  MRR/SDR  data,  the  probability  of 
corrosion  being  external  on  cacli  element  type  is  summarized 
in  Table  5. 

To  evaluate  this  function  during  the  simulation,  a 
random  number  from  a uniform  distribution  is  selected  for 
each  corrosion  occurrence.  If  t ho  random  number  is  between 
zero  and  the  probability  established  for  the  subject  ele- 
ment, the  corrosion  is  treated  as  being  external;  otherwise, 
it  is  treated  as  being  internal. 
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TABLE  5. 


PROBABILITY  OF  A CORROSION  DEFECT  BUI  NO  EXTERNAL 


To  t a 1 

Occurrences 


Fuselage 


Door  frame  13 
Floor  bourn  17 
Koel  beam  S 
Main  frame  22 
Stringer  257 
Pressure  web  0 
Window  frame  3 

Wing 

Access  frame  0 
Rib  3 
Spar  57 
Stringer  43 

Wing  Center  Section 

Rib  0 
Spanwise  beam  18 
Stringer  92 


* estimated 


lixte  rnu  I 

Occurrences 


8 

C) 

0 

2 

1 H 4 
0 
3 


n 

2 

id 

2 5 


0 

6 

32 


Probali  i 1 i i y of 
Corrosion  Being 

External  _____ 


0 . 0 1 5 
{).  352 
0.  1)01) 
0, 01)1 
n . 7 I (I 

0.  050* 

1 . non 


n.615* 

0 . 067 
[1.17  5 

tl  .531 


(I,  1 75* 
0.  333 
0.348 


Crack  defects  are  treated  in  the  same  manner  as  cor- 
rosion defects,  with  the  probability  of  a crack  beinji  ex- 
ternal, as  summarised  in  Table  b. 


TABLE  6.  PROBABILITY  OF  A CRACK  DEFECT  BE INC  EXTERNAL 


Fuselage 


Tot  a 1 

Occurrences 


External 

Occurrences 


l’rohnb  1 1 i t v of 
Crack  being 
Ex  to  mu  1 


Door  frame 

82 

JD 

0 . 2 3 2 

Floor  beam 

fill 

2 i 

1!  .2  41 

Keel  beam 

10 

0 . (ISO 

Main  frame 

735 

211 

O . 0 3D 

Pressure  web 

DO 

1 

0 . 11  1 I 

Window  frame 

50 

2 8 

0 . 8 (id 

Stringer 

•>50 

2 4 4 

0 , 32  5 

ing 

Access  frame 

1 1 5 

74 

0 . (.4  3 

Rili 

2 84 

30 

(1 . 1 lilt 

Spar 

(ill  7 

3(i  4 

0 . 54  (i 

S t r 1 ngcr 

51  a 

3D  7 

o , (i  ;i 

lng  Center  Section 

Rib 

i 

0 

[i , lino 

Stringer 

2 02 

132 

(1 . [i  5 3 

Spnnwl.se  beam 

1 18 

0 3 

0,  53  4 
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As  was  the  case  for  corrosion,  this  function  is  eval- 
uated by  selection  a random  number  from  a uniform  distribu- 
tion for  each  crack  occurrence.  Jf  the  random  number  is 
between  zero  and  the  probability  established  for  the  subject 
element,  the  crack  is  treated  as  being  external  ; otherwise, 
it  is  treated  as  being  internal.  However,  in  the  case  of 
cracks,  SMFIi  has  the  facility  for  cracks  which  initiate  in- 
ternally to  become  external.  In  general  when  a crack  be- 
comes external,  there  is  a higher  probability  of  it  being 
detected  since  the  lowest  level  external  inspection  is  nor- 
mally lower  than  the  lowest  level  internal  inspection.  Tito 
point  at  which  an  internal  crack  becomes  external  is  defined 
as  a percentage  of  the  critical  crack  length.  This  percent- 
age is  an  input  parameter. 

S . 5 Inspection  Covers  Defect  Area 

This  function  requires  two  program  inputs:  one  input 
to  specify  the  most  frequent  inspection  interval  that  covers 
the  internal  portion  of  the  subject  element,  and  a second 
input  to  specify  t lie  most  frequent  inspection  interval  that 
covers  the  external  portion  of  the  subject  element.  It  is 
assumed  that  once  the  most  . frequent  interval  is  specified, 
the  element  i.s  covered  during  all  higher-level  Inspections. 
The  simulation  is  set  up  to  work  with  a four- level  inspec- 
tion program.  The  levels  are  those  commonly  referred  to  by 
tlie  ATA  as  A,  Preflight;  R,  Service;  C,  Phase;  and  D,  Over- 
haul. The  frequency  of  each  level  is  also  an  input  param- 
eter, which  can  be  varied  to  reflect  the  inspection  program 
of  different  air  carriers. 

5 . <>  He  foe  t Pound 

5.6.1  Crack  Defect  Pound 


Tli c probability  of  detecting  a crack  depends  on 
the  defect  sice  ami  the  inspection  level.  The  basic  form  of 
the  probability  equation  Is  defined  by  Davidson  in  Reference 
IP.  The  equation  form  is  an  exponential  with  three  con- 
stants: Ci,  |A | , and  ?,0,  and  one  variable,  £ . The  proba- 

bility P(i')  for  each  inspection  level  is  defined  as  follows: 

0 for  e £ 

fl-oxpf-Hj  U-P.on  ) for  l ■>  tQ  (?7) 


where 


C,  ■ The  maximum  probability  of  detecting  a crack  at  u 
1 given  inspection  level.  If,  for  instance,  Ci 
" 0.9,  then  1 out  of  10  cracks  would  be  missed 
regardless  of  its  length.  Cj  approaches  1,0  for 


overhaul  ami  special  inspections,  hut  is  signifi- 
cantly less  than  1,0  for  lower  level  inspections. 
This  parameter  is  determined  from  MRK/SUR  data  by 
comparing  the  number  ol  detects  that  were  found  at 
a particular  inspection  Level  with  the  number  of 
defects  that  should  have  been  found,  but  were 
missed.  The  number  of  defects  that  were  missed 
was  determined  by  the  number  of  cracks  that  were 
found  on  higher* level  inspections  but  could  have 
been  detected  at  the  Lower  level  because  of  their 
size  and  location.  The  specific  defect  size  and 
location  for  each  inspect  ion *1 evel  was  determined 
from  the  Inspector's  Survey  and  from  the  inspec- 
tion level  definition.  The  criteria  for  each 
inspection  level  arc  as  follows: 

Proflight  - all  cracks  larger  than  0.551  inch 
and  located  on  the  lower,  exterior  surfaces  of 
the  wing  and  fuselage  could  have  been  detected  at 
th I s 1 eve  1 . 

Service  all  cracks  larger  than  0 . ■!  1 0 Inch 
and  located  on  the  exterior  of  the  upper  and  lower 
wing  surfaces,  the  lower  fuselage  surface,  and  the 
fuselage  near  the  wing  root  could  have  been  detec- 
ted fit  this  level. 

Phase  - all  cracks  larger  than  0.2(t(i  i licit  and  lo- 
cated on  the  wing  and  fuselage  exterior  surfaces, 
cargo  compartment , or  accessible  areas  of  the  wing 
interior  could  have  been  detected  fit  this  level. 

Overhaul  - all  cracks  larger  than  (i.LM  inch  in  all 
locations  could  have  been  detected  at  this  level. 

Because  of  the  very  few  MRK/SDR's  that  contained 
information  on  corrosion  size,  it  was  .assumed  that 
Cj.  for  corrosion  was  equal  to  for  cracks  fit  the 
same  1 eve  1 of  I nspoe t i on . 

B-  ■ The  shape  parameter  that  determines  how  quickly 
1 the  probability  of  detection  approaches  the  max- 
imum probability,  d..  ts . is  also  determined  from 
MRR/SDR  data  hy  plotting  the  cumulative  distri- 
bution of  defects  detected  at  a given  Inspection 
1 ovel  . 

£ » The  smallest  detect  that  will  he  detected.  This 

parameter  was  determined  from  the  data  in  Volume 
111,  the  Inspector's  Survey,  by  taking  the  loga- 
rithm of  the?  responses  at  a given  inspection  level 
and  computing  the  mean  value.  The  mean  of  the 
logarithms  was  used  because  a plot  of  the  re- 
sponses indicated  that  they  closely  followed  a 
log  no  rtitfi  1 d i s t r i but  ion. 

1)5 


Having  determined  the  probability  of  detection  for  a given 
defect  size,  the  determination  of  whether  or  not  a defect  was 

found  is  based  on  Monte  Carlo  methods.  A random  number  is  gen- 
erated from  a uniform  distribution  and  if  that  number  is 
less  than  or  equal  to  the  probability  of  detection,  the 
defect  is  considered  to  he  detected;  otherwise,  it  is  con- 
sidered to  be ' undetected . 

If  the  defect  Is  detected,  the  element  Is  repaired  or 
modified  depending  on  previous  simulation  decisions.  If  the 
defect  remains  undetected,  the  simulation  continues  to 
conduct  inspections  until  the  defect  is  detected,  a struc- 
tural failure  occurs,  or  the  aircraft  is  retired  from  the 
fleet . 


5,6,2  Corrosion  Defect  found 


The  probability  of  detecting  a corrosion  defect  is 
determined  similarly  as  the  probability  of  detecting  a crack 
defect.  The  basic  equation  is  the  same,  except  that  the 
corrosion  area,  a,  is  substituted  for  the  crack  length,  , 
and  the  minimum  detectable  area,  a0,  Is  substituted  for  the 
minimum  crack  length,  t0.  The  equations  for  the  probability 
of  detection  of  corrosion  then  become  the  following: 

Pre flight:  A Level 

Pin)  » 0.27M1.0  - exp  f - 0 . 51  2 (a  ••  0 . '.)'.»(»)  ] i (421 


Service:  B Level 

P (a  ) = <l.4n2(l.l)  - exp  | - 0 . 5 I 3 ( a - U . 63  51  | > 


Phase:  C Level 

P ( a ) =>  O.lvI.Kl.O  - oxpl  • (1.  (i8(.(a  - ().5(i,V)|l  (44) 

Overhaul  : I)  Level 

P (a ) * 0.  DIMM  1.0  - exp  | - 0 . 72  5 (a  - 0.351)]}  (45) 


These  equations  are  Illustrated  in  Figure  4o. 

Again,  a Monte  Carlo  method,  using  a uniformly 
distributed  random  number,  is  employed  to  determine  whether 
a defect  is  actually  found  or  not. 


6 7 


OirJvMf  ifli  tfiMi  i fcinifafcBt  


‘“ton  mi1 


1,1  VHI.^JNSIMUITfON  _ jiVUIOIAIIJ 


I4 


I.IVI-I.  INSI’liCTlON 


VI! I.  INSI'IUITlfl 


A"  l.liVlil,  INK 


OMkOSKtN  AIWA  (In,) 

l-‘  1 >>  11  rt'  lo.  Probability  of  Corrosion  Detection  During  a 
Periodic  inspection 


5 • 7 Inspect  Ion  frequency  of  li lenient  Decreased 

Inspection  Intervals  are  normally  extended  in  structural 
areas  where  few  defects  have  boon  found.  As  the  overhaul  and 
phase  Inspections  are  conducted  on  each  of  the  ten  high-tlme 
aircraft  in  the  fleet,  the  time  of  detection  and  the  number 
0 1,  delects  detected  are  recorded.  If  no  defects  are  found  on 
any  of  the  fen  high-time  aircraft  during  one  D-levol  interval, 
then  the  overhaul  Inspection  interval  is  extended.  The  amount 
ol  decrease  depends  on  the  particular  element  and  is,  there- 
fore, an  input  parameter. 

Inspection  inteival  extensions  apply  only  to  the  phase 
and  overhaul  inspections. 

Time  to  Next  Inspection  l.ess  Than  Mine  to  flemont 

Failure 

Tills  function  compares  the  time  to  structural  failure 
with  tho  time  of  the  next  periodic  inspection.  if  the  time 
to  the  next  inspection  occurs  first  and  the  aircraft  is  not 


scheduled  for  retirement,  then  the  inspection  routine  is 
repeated.  It'  the  time  to  failure  occurs  first,  the  logic 
continues  to  Block  5.10, 

5 . 9 Aircraft  Service  Life  hxpircs  Pr  ior  to  Next. Inspect,  ion 

This  function  simply  compares  the  service  life  of  the 
aircraft,  as  stated  hv  the  manufacturer  or  as  determined  by 
the  user  from  service  experience,  with  the  number  of  flight 
hours  that  the  aircraft  has  accumulated  and  the  number  of 
flight  hours  until  the  next  periodic  Inspection  is  sched- 
uled, If  the  accumulated  flight:  hours’  plus  the  hours  to  the 
next  Inspection  are  equal  to  or  greater  than  the  aircraft 
service  life,  then  the  aircraft  is  retired  from  the  fleet. 
If,  however,  the  service  life  has  not  expired,  the  logic 
returns  to  Block  5.1  and  the  inspection  routine  is  repeated 
until  a defect  is  found,  a structural  failure  occurs,  or  the 
aircraft  is  retired  from  service. 

5.10  De  l e t e A i r c r a f t 1;  r um  i Jo  e t 

Aircraft  are  deleted  from  the  fleet  as  their  service 
life  expires  or  as  they  experience  structural  failures,  The 
aircraft  service  life  is  an  input  parameter  that  can  he 
varied  to  determine  the  effect  on  safety  of  extending  the 
service  life  with  or  without  improving  the  design  require- 
ments . 

5.11  Last.  Aire  raj  t in  fleet. 

This  function  sirply  keeps  track  of  the  number  of 
aircraft  deleted  from  t he  fleet  and  compares  that  number 
with  the  number  of  aircraft  originally  produced.  When  the 
last  aircraft  lias  been  deleted,  SAII:K  returns  to  Block  1 . U 
and  repeats  the  simulation  with  the  next  element. 


(i . Dot  a ih'il  Description  of  hi  och  h.  0 , Repair  I:  lenient  to 
0 r i g 1 it  a 1 Strength  f see  i;  inure  4 7) 


figure  4 7.  Detailed  Plow  Diagram  - Repair  lilement  to 
Original  Strength 


0,1  Modification  Pend i ng 

As  dismissed  in  Section  S . 3 , if  there  is  a structural 
modification  pending  on  a particular  aircraft,  the  SATPii 
Logic  represents  the  installation  of  the  modification  as 
occurring  during  either  an  overhaul  inspection  or  the  repair 
of  an  in-service  defect.  If  there  are  future  in-service 
defects  projected  for  the  aircraft  at  the  time  of  installa- 
tion, these  projections  are  cancelled  and  new  ones  are 
generated  from  ’their  corresponding  distributions  as  though 
the  aircraft  were  just  entering  service.  If  the  inspection 
intervals  had  beer.  pr>'\  leusly  decreased,  t hey  are  returned 
to  tbeir  initial  values  as  each  aircraft  is  modified. 

(> . 2 Ac  tun  1 Pa  t i gue  1.  i t'e  o f Mod  i fled  P t ement 

When  a structural  modification  is  installed  on  an 
aircraft,  the  actual  fatigue  life  of  the  modified  element  is 
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determined  in  the  same  manner  as  when  the  aircraft  entered 
service.  The  actual  life  is  drawn  from  a stat  ist  ical  dis- 
tribution of  fatigue  lives  about  the'  actual  average  fatigue 
life  of  the  element  type.  If  the  modification  has  not  been 
fatigue  tested,  the  actual  average  fatigue  life  of  the 
element  type  is  statistically  determined  as  it  was  at  the 
start,  of  the  simulation  in  Block  1.4.  If  the  modification 
has  been  fatigue  tested,  the  actual  average  fatigue  life  of 
the  element  type  is  assumed  to  be  the  same  as  that  predicted 
by  analysis. 

6. 3 Rep  a i red to  Original  Actual  Average  Juit  igue  1,1  fe 

When  it  has  been  determined  that  a repair  is  required 
and  if  there  is  no  modification  pending  on  the  element,  the 
element  is  repaired  before  any  more  flight  hours  arc  allowed 
to  accumulate  on  the  aircraft.  All  defects  present  at  the 
time  of  repair,  whether  detected  or  not  during  the  scheduled 
inspection,  arc  assumed  to  he  repaired  at  this  time.  It  is 
also  assumed  that  the  element  strength  is  restored  in  its 
original  static  strength,  However,  in-service  defects 
(fatigue  cracks,  service  damage,  and  corrosion)  predicted  to 
occur  after  the  repair  is  accomplished  are  not  affected  by 
the  repair  process  and  they  are  allowed  to  occur  at  their 
originally  detenu  ined  times. 

0 . 4 Actuu_l_  fat  i gue  hi  fe  of  Repa.i  r li  1 e men  t 

As  discussed  in  Section  t> , 3 , previously  projected 
fatigue  cracks  that  have  not  been  initiated  at  tiie  time  of 
repair  are  unaffected  and  retain  their  original  initiation 
times.  Those  cracks  that  arc  repaired  have  new  times  to 
crack  initiation  determined  in  tiie  same  manner  as  when  the 
aircraft  entered  service;  that  is,  from  a fatigue  life 
distribution  reflecting  basic  fatigue  scatter  and  load 
environment  variation,  times  are  randomly  drawn  about  the 
element  average  fatigue  life. 


Detailed  Description  of  Block  7.0,  Special Inspec t ion 

and  I n c r o7i I n spec t i on  _!■' r e<|uc ncy  (see  Figure  48j 


Figure  48. 


Detailed  Flow  Diagram  - Special  Inspection  and 
Increase  Inspection  Frequency 
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7.1  [ ns  poet  i on  I ^requeue  y of  Lloment  Dos  i gn  I ncreasod 

The  frequency  at  which  certain  inspections  are  con- 
ducted is  increased  when  it  has  been  determined  that  the 
present  frequency  is  not  adequate.  The  percentage  increase 
in  frequency  is  an  input  parameter  that  depends  on  the 
safety  criticality  of  the  element.  Three  criteria  arc  used 
to  determine  whether  a frequency  increase  is  necessary.  If 
any  one  of  the  three  are  satisfied,  the  frequency  is  in- 
creased, The  throe  criteria  arc  as  follows: 


(.  1 1 A crack  greater  than  fail-safe  length  is  detected. 

(2)  Large  cracks  are  detected  in  an  .individual  ele- 
ment, such  that  the  sum  of  the  lengths  of  the 
cracks  present  plus  the  projected  growth  of  the 
largest  crack  through  the  next  inspection  interval 
results  in  a one  half  reduction  in  the  I’m  i I safe 
st  rong  t li . 

(?)  Small  cracks  are  detected  in  the  same  element  on 
numerous  aircraft,  such  that  the  total  strength 
reduction  resulting  from  all  of  the  cracks  divided 
by  the  number  of  aircraft  in  the  fleet  equals  2 1) 
percent  of  the  original  fail-safe  strength  of  an 
i nd  ividua 1 e 1 cment . 


The  first  two  criteria  deal  with  defects  that  very 
seldom  occur,  but  when  they  do  there  is  a high  probability 
of  an  aircraft  accident  or  of  extensive  unschedu I ed  main- 
t o nance . 


The  third  criterion  deals  with  the  potential  safety 
hazard  resulting  from  the  very  small  but  finite  possibility 
of  the  occurrence  of  a flight  load  that  exceeds  the  design 
strength  along,  with  the  greatly  increased  probability  that 
the  strength  of  the  element  Ikis  been  slightly  reduced  be- 
cause o f a sin  a I I c rnc  k . 


7 . 2 Sc  bed  ti.1  el  mined  i at  e Sp  e c i a 1 1/  I eet  I aspect  i on 

The  criteria  defined  in  Section  ’".I  are  also  used  to 
determine  whether  a special  fleet -wide  inspect  ion  is  re 
quirod.  When  a special  inspection  is  called  for,  the  sub- 
ject element  is  carefully  inspected  on  every  aircraft  in  the 
fleet.  If  a special  inspection  is  not  called  for,  the  logic 
continues  to  Block  N . U wiiere  it  is  determined  whether  or  not 
a modification  to  the  element  is  required. 
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7.3  lie feet  Present 

This  function  is  ;i  simple  yes-or-no  indicator  that 
determines  whether  a crack  defect  has  initiated  on  the 
element  being  covered  by  the  special  inspection.  If  no 
cracks  are  present,  the  logic  continues  immediately  to 
mock  8.0. 

7.4  Defect  found 


During  a special  inspection,  the,  probability  of  finding 
a defect  is  significantly  improved  over  the  probability  of 
finding  a defect  during  a regularly  scheduled  inspection. 
This  is  duo  to  the  fact  that  the  locution  and  nature  of  the 
defect  are  reasonably  well  specified  before  the  inspection 
is  conducted.  The  probability  of  detecting  a crack  during  a 
special  inspection  is  determined  from  the  equation 


P(l)  ■ 0.090(1.0  - exp [-0.971  ( t - 0.102)]}  (46) 

And  the  probability  of  detecting  corrosion  is 

P(a)  a 0.9991 1.0  - exp  l -0.971  (11  - 0.180)]}  (47) 


where 


l » length  of  the  crack  present 
a * area  of  corrosion  present 

As  was  the  case  for  the  probability  of  detection  curves  in 
Section  5.6,  these  equations  were  determined  from  the  de- 
fects reported  in  the  MRR/SDR  data  and  the  Survey  of  Inspec- 
tors; these  equations  are  illustrated  in  figures  49  and  50. 
The  MRR/SDR  data  are  presented  in  Volume  III  of  this  report. 
The  defect  is  considered  to  be  found  when  a random  number 
generated  by  the  simulation  is  equal  to  or  less  than  the 
probability  generated  from  the  above  equations  for  cracks  or 
corrosion . 

7.5  Inspection  frequency  Increased  because  of  Special 
1 1 \~s l •> e c t To n R o 'suit's  ' - - — 

After  a special  fleet  Inspection  is  completed,  the 
magnitude  of  the  defects  found  is  compared  with  the  first 
and  second  criteria  defined  in  Section  7.1,  and  a second 
reduction  in  inspection  frequency  may  be  instituted  as  a 
result  of  the  special  inspection,  for  the  inspection  in- 
terval required  by  the  second  criterion  in  Section  7.1,  the 
interval  most  recently  set  by  flock  7.1  is  used. 
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PRESENT  SERVICE  EXPERIENCE 
WARRANTS  MODIFICATION 


DEFECTS 


PRESENT  SERVICE  EXPERIENCE 
WARRANTS  ADDITIONAL  MODIFICATION 


MODIFICATION 
TESTED  (INPUT) 


PREDICTED  AVERAGE 
MODIFIED  FATIGUE 
1. 1 PE 


ACTUAL  AVERAGE 
MODIFIED  FATIGUE 
LIFE 


INSTALL 


ACTUAL  AVERAGE  MODIFIED 
FATIGUE  LIFE  SET  EQUAL 
TO  PREDICTED  LIFE 


Figure  SI.  Detailed  Flaw  D Ingram  - Develop  Modification 
Because  of  Service  lixperience 
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8 . 1 11 1 ement  P rev  i ous  l)'  Mod  i f i od 

As  a simple  variable,  this  I unction  .indicates  whether 
an  element  was  previously  modified.  If  the  element  was  not 
modified,  Block  8.2  examines  present  service  experience  and 
makes  a decision  on  the  development  of  a first  modification. 
If  the  element  was  previously  modified,  Block  8.5  examines 
present  service  experience  and  makes  a decision  on  the 
development  of  an  additional  modification. 

8.2  Present  Service  Experience  Warrants  beve  l opine n t'  of 

a ModTFlcatMbn 

The  S A l 111  logic  bases  the  decision  on  whether  or  not  to 
develop  a structural  modification  because  o I'  service  exper- 
ienee  solely  on  economic  considerations.  Since  the  economic 
parameters  considered  depend  on  element  type  and  are  subject 
to  change  with  time,  they  are  necessarily  part  of  the  input 
data.  Values  for  most  of  these  parameters  are  related  to 
the  inspection  level  at  which  the  maintenance  action  is 
performed.  This  relationship  is  required  because  the  com- 
plexity of  the  elements  generally  increases  with  higher 
inspection  levels.  The  economic  parameters  identified  and 
the  values  currently  used  are  listed  in  Table  7.  The  valuers 
shown  in  Table  7 were  determined  from  a report  t Reference 
20)  given  at  the  /VIA  Maintenance  Conference  held  September 
20  to  28,  l!)(nS,  The  estimated  values  are  based  on  the 
relative  complexity  of  each  inspection  level. 


TABl.li  7.  IldONliMIO  PARAMIITIIRS  Rliqil  I lUll)  PUR  1 ) II V I i 1,0 1* I Nti  A 
SPRY  I Cli  MODI  P I OAT  I ON 
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Man  -Mrs/  1 asp . 
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The  decision  to  develop  a modification  is  made  by 
comparing  the  cost  per  flight  hour  of  the  modification  with 
the  repair  cost,  per  flight  plus  the  increased  inspection 
cost  per  flight  hour.  The  modification  cost  per  flight  hour 
is  found  by  dividing  the  total  fleet  modification  cost  by 
the  remaining  service  life  of  the  fleet.  The  repair  cost 
per  flight  hour  is  found  by  dividing  the  total  fleet  repair 
costs  since  the  Last  modification  by  the  fleet  flight  time 
since  the  last  modification.  The  increased  inspection  cost 
per  flight  hour  is  found  by  dividing  t lie  projected  increased 
inspection  costs  by  the  remaining  service  life  of  the  fleet. 
A modification  is  justified  when 


111;  PAIR 


TNSP 


where 


'Rl;  PAIR 


= modification  cost  per  flight  hour 
» repair  cost  per  flight  hour 


! 1 ns p “ increased  inspection  cost,  per  flight  hour 


8.3  Pres e n t.  3 e r v 1 c e _ !; xpor i enco  Warrants  Development  of 

AcfdT t T o hal " ModTTTc a t urn 


The  decision  to  develop  an  additional  modification 
because  of  service  experience  is  based  on  the  same  economic 
parameters  considered  in  Block  8.2.  Currently  the  economic 
parameter  values  used  are  the  sumo  as  those  for  the  first 
modification  decision;  those  values  are  listed  in  Table  7. 

8 . 4 Mod  i locat  i on  Tested 

This  function  indicates  whether  or  not  a modification 
has  been  tested  before  its  incorporation  into  the  fleet.  It 
is  an  input  parameter  that  is  determined  for  each  element. 

8 • 5 Predicted  Average  Modified  Fatigue  Life 

if  a decision  has  been  made  to  develop  a modification, 
it  is  assumed  that  the  mod  i f i cat  i on  .is  again  designed  to  the 
predicted  average  fatigue  life  of  the  original  design  as 
defined  in  Block  1.3. 

8 . 6 Actual  Average  Modi  fled  Fatigue  Life 

If  the  modification  is  not  fatigue  tested,  it  is  sub- 
ject to  the  same  type  of  variation  between  uctual  and  pre- 
dicted fatigue  life  as  was  the  original  design.  To  deter- 
mine the  actual  average  fatigue  life  of  the  modification,  a 
random  draw  is  once  more  made  from  the  log-normally  dis- 
tributed correction  factor  described  in  Section  1.5.  Al- 
though the  form  of  the  correction  factor  distribution  is  the 


same  as  it  was  for  the  original  design,  it  is  assumed  that 
the  accuracy  of  the  modification  design  is  improved.  This 
increased  accuracy  is  accounted  for  by  decreasing  the  stan- 
dard deviation  and  increasing  the  mean  in  the  above  disti  i 
but  ion.  The  standard  deviation  is  reduced  by  15 « and  the 
mean  is  increased  by  the  quantity  0 . 1 S (1.0 -mean). 

8.7  Actual  Average  Modified  fatigue  l.il'e  Set  Hqual  to 

FTiUTcted'Tlte' 

If  the  decision  was  made  to  fat.ig.ue  test  the  modifi- 
cation, it  is  assumed  that  the  actual  average  fatigue  life 
of  the  modification  will  attain  its  predicted  life  or  he 
redesigned  and  retested  until  it  does. 


III.  SUMMARY 

The  implementation  of  the  SAII'li  logic,  as  described  in 
Section  11,  retiuircd  incorporating  decision  functions  to 
realistically  simulate  the  evaluation  process  which  es- 
tablishes or  modifies  the  structural  inspection  intervals 
for  commercial  jet  transport  aircraft.  Some  of  these  de- 
cision functions  were  based  on  previous  research  as  refer- 
enced in  this  volume,  and  others  were  developed  during  the 
current  study.  A partial  list  of  the  lHttcr  include  the 
following: 

(1)  Probability  of  a production  defect  occurring. 

(2)  Service  damage  occurrence  rates. 

(3)  Corrosion  occurrence  rates. 

(4)  Probability  of  corrosion  occurring  in  a stress 
concentration. 

(5)  Probability  of  crack  detection  during  a periodic 
or  a special  inspection. 

(6)  Probability  of  corrosion  detection  during  a per- 
iodic or  a special  inspection. 

(7)  Pressurization  load  exceedances  per  flight  hour. 

The  results  of  the  demonstration  computer  runs  indicate 
that  the  implemented  SAIFK  logic  can  successfully  simulate 
and  quantify  the  evaluation  process  for  the  establishment  or 
revision  of  the  structural  inspection  intervals  for  commer- 
cial jet  transport  aircraft. 
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Figure  52.  Flow  Diagram  Showing  Major  Aspects  of  SA I F Li  Logic 
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I Figure  54.  Detailed  Flow  Diagram  - Develop  Modification 

I Because  of  Fatigue  Test  Failure 
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Figure  SS.  Detailed  Flow  Diagram  - Reduce  Fatigue  Life 

Because  of  Production,  Service,  or  Corrosion  Defects 


Figure  56.  Detailed  Flow  Diagram  * Reduce  Strength  Because 

of  Crack/Corrosion  Growth/Predict  Time  to  Failure 
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Figure  57.  Detailed  Flow  Diagram  - For  iodic  Inspection 
of  Iilements/Aircraft  Deleted  from  Fleet 
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Figure  60.  Detailed  Flow  Diagram  - Develop  Modi  Tit;  at  ion 
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